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CHARTS AND APPROXIMATE FORMULAS FOR THE ESTIMATION OF AEROELASTIC EFFECTS 

ON THE LOADING OF SWEPT AND UNSWEPT WINGS ' 


By Franklin W. Diedbrich and Kenneth A. Foss 


SUMMARY 

Chartu and approximate formulas are presented for the 
eshmation of aeroelastic effects on. the spanwise lift distri- 
bution, llft-curre slope, aerodynamic center, and damping 
in roll of swept and unswept wings at subsonic and supersonic 
speeds. Some design considerations brought out by the results 
of this report are discussed. 

INTRODUCTION 

A Unowloclgo of Ihc spanwise lift distribution and of some 
of tlie aei'odynaniic pai-anielers associated with it is recpiired 
for flic desif>:n of a wing structure. Under certain conditions, 
sucli as liigli dynamic i)ressures, thin wings, swept wings, 
or wings designed for low wing loadings, the spanwise lift 
disti'ihul ion may he affected to a significant e.Ntent by 
aeroelastic effects, because a wing which carries a certain lift 
necessarily deforms under that lift. If the angles of attack 
along the span are changed as a result of this deformation, the 
lift cai'ried by the wing is changed as well; in turn, this 
cliangi' in lift causes a change in the deformation of the wing 
and hence another change in lift, and so on, until an equilib- 
rium condition is reached. The changes in the. magnitude 
and the distribution of the lift are reflected in changes of the 
wing lift-curve slope, the wing bonding and rolling moments, 
tlie spanwise center of pressure of the lift, and, on a swept 
wing, llic longitudinal center of pressure. 

Inasmuch as the lift produced by a given change in angle 
of attack is pi-oportional to the riAmamic pi-essure, the various 
aeroelastic effects tend to increase with dynamic pressure. 
In fact, for certain wings a sufficiently large dynamic. pressure 
may produce a condition of instability in which the change in 
lift caused by deformation is greater tlian the amount of 
lift recpiired to produce the deformation, so that a given 
defoi-ination will tend to increase until the structure fails. 
This phenomenon is aeroelastic divergence; since it involves 
only torsional defoi-rnations in the case of unswept wings, it is 
often referi-ed to as torsional divergence. 

Several methods are available for calculating these efl'ccts 
(ref. 1, for instance), but since these effects depend on the 
structiwal characteristics of the wing, which are not ac- 


curately known in advance of its design, the relatively large 
amount of time required for even the most efficient of these 
methods militates against their use in connection with 
preliminary design calculations. A need e.vists, therefore, 
for means of estimating some of the more important aero- 
elastic effects on the spanwise lift distribution quickly and 
with an accuracy that is sufficient for preliminary design 
purposes. 

Charts and appro.ximate formulas are presented in this 
report for estimating the changes in spanwise lift distribution, 
lift-curve slo])e, wing rolling-moment coefficient, spanwise 
center of pi-essure, and aerodynamic center occasioned by 
aeroelastic action of swept and unswept wings at subsonic 
and supersonic speeds. Also included arc summary charts 
which indicate whether the various aeroelastic phenomena 
considered are likely to affect any given design. By means 
of these charts the conventional procedure of designing a 
wing on the basis of certain strength criteria, checking it for 
aeroelastic phenomena, and then reinforcing it. when neces- 
sary, to meet the stiffness requirements imposed by these 
phenomena can often be simplified greatly, inasmuch as 
the effect of some of these phenomena can be estimated in 
advance of design. 

The use of the charts is described in the section headed 
“Calculation of the ^'arious .Aeroelastic Phenomena,” and 
some considerations involved in the selection of the aero- 
dynamic, structural, and geometric parameters are discussed 
in some detail in the section headed “Selection of Param- 
eters.” These two sections, as well as tJic sections headed 
“Illustrative E.xample” and “Preliminary Survey of Aero- 
elastic Behavior,” are likely to prove of greatest interest 
at a first reading of this report. The various parts of the 
section headed “Discussion” are concerned with the limita- 
tions of the charts, with the light they shed on such practical 
design problems as the relative significance of strength and 
stiffness as design criteria, with efficient ways of stiffening 
a. wing that is strong but not stiff’ enough, and with the 
achievement of aeroisoclinic conditions. 

A brief description of the calculations (based on refs. 1 
and 2) used in preparing the charts is contained in the 
appendixes. 


* Urcviously relensod ns XACA 1'N 2G05i, “Clinrts and Approviinalc Kormiilas for tho Estimation of Aoroclastic EiTecis on tin* Loading of Swept and Unswepi Wings,” by Franklin W. 
Diodorieh and Kunnclli .A. Fos-?. lOW, 
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SYMBOLS 

aspect ratio, b'-fS 
swepl-span aspect ratio, A/cos-\ 
cross-seetional area of the (assumed) single 
torsion cell, sq in. 

distance from leading edge to section aerody- 
namic center, fraction of chord 
distance from leading edge of mean aei'odynainic 
chord to wing aerod^mamic center, fraction of 
mean aerodynamic chord 
wing span, in. 

wing span less width of fuselage, b — w, in. 
wing-root bending-moment coefficient, 4Afr/qSb 
lift coefficient of wings alone (e.xclusive of fuse- 
lage), LJqS 

wing lift-curve slope per radian 
efi’ective lift-curve slope per radian (defined in 
eqs. 2 and 4) 

rolling-moment coefficient on both wings alone 
(exclusive of fuselage). Rolling momonl/qSb 
wing-root twisting-moment coefficient, 2Tr/qScr 
chord (measured perpendicular to elastic axis), 
in. 

average chord, in. 

section lift-curve slope per radian 
mean aerodynamic chord (parallel to plane of 
symmetry), in. 

Young’s modulus of elasticity, Ib/sq in. 
distance from leading edge to elastic axis, frac- 
tion of chord 

dimensionless moment arm of the section lift 
about the elastic a.xis, e — a 
effective or average dimensionless moment arm 
allowable bending stress, Ib/sq in. 
root-stiffness function 
allowable shear stress. Ib/sq in. 
structural weight function 

dimensionless parameters used in approximate 
formulas for angle of attack due to aero- 
elastic deformation 

dimensionless functions of the distance along 
the span used in approximate formulas for 
angle of attack due to aeroelastic deformation 
modulus of rigidity, Ib/sq in. 
wing thickness, in. 

section bending moment of inertia, in.'' 
section moment of inertia in torsion, in.^ 
dimensionless parameters used in approximate 
formulas for dimensionless dynamic pres- 
sures at divergence 

dimensionless sweep parameter, ~ tan A 

c,c, {EI)r 

lift of both wings alone (exclusive of fuselage), lb 
lift per unit distance along span, Ib/in. 


bending moment about an axis porpendieular 
to elastic axis, in-lb 
free-stream Mach number 
design load factor 
rolling angular velocity, radians/sec 
dynamic pressure, Ib/sq ft 
dimensionless dynamic pi’cssuic, 
q CL„fiCr-sr cos A 
l44 WJ)r 

dimensionless dynamic jiressure, 
q sin A 

144 (E7b 
total wing area, sq in. 

distance along elastic axis measured from wing 
root, in. 

dimensionless distance along elastic axis, •'f/.v, 
distance from root to center of pressure of lift 
along elastic axis, in. 

dimensionless distance from root to center of 
])ressure of lift, sjs, 

accumulated torque about elastic axis, in-lb 
thickness of most highly stressed element of 
skin, in. 

thickness of equivalent skin which includes the 
material in stringers and spar flanges, in. 
distributed torque due to inertia loading, in-lb/in. 
free-stream velocity, ft/sec 
design gross weight of airplane, lb 
weight of primary structure of both wings, lb 
weight of both wings exclusive of fuselage, lb 
width of fuselage, in. 

weight of primaiy load-cariying structure per 
unit distance along span, Ib/in. 
lateral coordinate, in. 

7/ 

dimensionless lateral coordinate, ^ 

lateral distance to center of pressure 
angle of attack in a plane parallel to plane of 
symmetry, radians 

angle of local dihedral, radians; or spanwise 
slope of normal displacement of elastic axis 
density of the material of the primaiy structure 
(or an equivalent density in the case, of sand- 
wich construction), Ib/cu in. 
lateral distance measured from wing root, 
w . 

y~2’ 111 . 

dimensionless lateral distance, ^7^;^ 

factors defined in table 1 

factors defined in equations (15a) and (lob) 

ratio of lift-curve slopes, 

angle of sweepback at elastic axis 
wing taper ratio, cjcr 
free-air density, slugs/cu ft 
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ip angle of slnictiiral twist in planes peipendicular 

to clastic axis, radians 
cj ti]) stiffness ratio, (EI),l(EI)i^^ 

H,v,T dimensionless parameters used in apiiroximate 

formulas for lift, root bending moment, and 
root twisting moment 

Subscripts; 

cx constant stress 

I) at divergence 

c effective 

g geometric (due to airplane attitude or built-in 

twist) 

i inertia 

0 rigid wing (for 5 = 0) 

/• at wing root 

.s- structural (due to structural or aeroelaslic 

deformation) 

1 at wing lip 
Superscripts: 

M due to bending moment 

T due to torque 

r due to root bonding 

ip due to root twist 

USE OF THE CHARTS AND APPROXIMATE FORMULAS 

SUMMARY OF METHOD AND SCOPE OF CALCULATIONS ON WHICH THE 
CHARTS AND APPROXIMATE FORMULAS ARE BASED 

Altliough a detailed understanding of the method and 
sco]ie of the calculations on which the charts of this report are 
liascd is not essential to the use of the charts, a brief account 
of these matters is given, primaril}’" to aid in the appreciation 
of the limitations of the charts. The method is described 
more fully in appendix A. 

Most of the calculations on which the charts are based 
were made by the method of refei'cnce 1, which consists in 
solving the difi’ei-ential equations descriptive of an elastically 
deformed wing under aerodjniamic loading by numerical 
methods employing matrix techniques. Treated by this 
method were wings with four taper ratios X (1.0, 0.5, 0.2, and 
0), two types of stiff'ne.ss distributions (one propoilional to 
the fourth power of the chord and one dictated by constant- 
stress considerations), and four values of a sweep parameter 

A- at sevei-al values of the dynamic-])ressure ratio Cal- 
culated for each case were the dynamic pressure at divei'gence 
and the changes due to aeroelastic action in spanwiso lift 
distribution, total wing lift, root bending moment, rolling 
moment, and spanwise center of pre.ssure of the lift. For 
the wings of constant chord and constant stiffness, calcula- 
tions wore also ])erfoi-med for six values of k by a method 
which is an extension of that of reference 2 and consists in 
solving the tliff'erenf iai ecpiations exactly for these relatively 
simple cases. 

Some ai)])roximations have been made in the calculations 
concerning the aerodynamic induction effects, the root rota- 
tions, and the stiffness distributions, primarily in order to 
bold the number of variables considered in the analy.sis to a 


minimum and to make the results more generally applicable. 

Aerodynamic induction effects at subsonic speeds are 
taken into account by an overall reduction of the strip- 
theory loading and, in the matrix calculations, by rounding 
off the strip-theory loading at the tip (see refs. 1 and 2); for 
supersonic speeds, strip theory is used with a small reduction 
at the tip in the matrix calculations. This approximation 
has made it unnecessaiy to consider explicitly the effects of 
aspect ratio, sweep, and Mach number on the rigid-wing lift 
distribution; the effects of these parameters on the total lift 
and on the aeroelastic increment to the lift distribution have 
been taken into account. 

The rigid-body rotations imparted to a swept wing by its 
triangular root portion vaiy among different designs in a 
largely unpredictable manner. They have therefore been 
taken into account only by the use of an effective root, the 
selection of which in any given case is discussed briefly in a 
subsequent section. 

The spanwise distributions of the bending and torsional 
stiffnesses depend on the detailed design of the wing and 
cannot be generalized easily. The stiffness distributions 
used in the calculations of aeroelastic effects were obtained 
from the constant-stress concept outlined in appendix B, 
which constitutes an effort to relate the stiffness of a wing to 
its strength on the basis of the following assumptions: 

(1) The level of combined bending and torsional stresses 
is constant along the s]ian. 

(2) The structure is designed for combined bending and 
torsional stresses in such a manner that the sum of the ratio 
of the actual to the allowable bending stress and the ratio of 
the actual to the allowable torsion stress is equal to unity 
when the margin of safety is zero. 

(:i) The structure is of the thin-skin, stringer-reinforced 
shell type and its main features do not vaiy along the span; 
for instance, the number of spars and their chordwise loca- 
tions are constant along the span. 

(4) At the design condition the spanwise distribution of 
the applied loading is proportional to the chord. 

Also used in the calculations were stiffness distributions 
which vary as the fourth power of the chord, as do those of 
solid wings and wings with geometrically similar cross sec- 
tions; as pointed out in a subsequent section, the results of 
these calculations can be used to estimate aeroelastic phe- 
nomena of some wings which have large cutouts or which 
for some other reason do not have stiffness distributions 
represented fairly closely by those of the constant-stress 
type. 

All caleulatious are based on the assumptions that twisting 
is resisted primarily by the torsion cells of the wing structure 
and that the wing deformations can be estimated by means 
of the elementary theories of bending and torsion about an 
elastic axis. 

SELECTION OF PARAMETERS 

Geometric parameters. — -The geometric parameters used 
in the analysis are defined in figure 1. The location of the 
effective root indicated in figure I is discussed in the section 
concerned with the structural parameters. 
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Aerodynamic parameters. — Tlic acrodyiiamu- parameters 
which enter the analysis are the wing lift-curve slope and 
the location of the aerodynamic center. Two lift-curve 
slopes are used at subsonic speeds: The wing lift-curve 
slope is used only in conjunction with additional lift 
distributions; for all other lift distributions — that is, those 
due to built-in twist, due to roll, or due to aeroelaslic twist — 
an efl'eetive lift-curve slope is used. Approximate 

values of lliesc parameters arc given for suberitieal speeds 
by the relations 




.1 cos .V 

c. 


vl -r 2 cos A 

'2tt 




cos A 

.1 + 4 cos A 
Jtt 


( 1 ) 


( 2 ) 


wliere C;^ is the lift-curve slo])e of the section iierpendieular 
to the quarter-chord line at a Mach number equal to Mo cos A. 
An appro.ximate value is given by 


Ci 


X 1 — COs’A 


(3) 


Equation (1) is given in reference .3 and shown to be appli- 
cable both to incompressible and to suberitieal compre.ssible 
How. Equation (2) is given in references 1 and 2 but without 
the term in the denominator. This term is intro- 

duced into equation (2) in order to c.xtend its aj)i)licability 


to compressible Hows in the same manner as that em|)loyed 
for the coefficient of damping in roll in reference 3. (The 
lifting-surface corrections given in ref. 3 for the wing lift- 
curve slope and foi' the coefficient of damping in roll Jtre Jtot 
inelufled in e<is. (1) and (2) because they are imiiortant 
primarily foi’ wings of very low aspect ratio, to which the 
method of the present report is not applicable.) 

At sujiersonic speeds (more specifically, for supersonic 
leading and trailing edges) both lift-curve slopes are approx- 
imately equal to the effective section lift-curve slope; that is 


where 




_ 4 cos A 

‘'"^~x-U„-cos- A— 1 


( 4 ) 


The ratio of the lift-eurve slojies and C[,^ is defined by 


so that for supersonic speeds k is equal to 1. 

The local aei'odynamic centers^are^assumed to .be a t^a 
constant fraedion of the chord from the leading edge, so that 
their distances from the leading edge (as fractions of the 
local chords) ajc all equal to the distance of the wing aero- 
dynamic center fj-oin the leatling edge of the mean aerody- 
namic cliord (as a fraction of the mean aerodynamic chord). 
'I'he monuMit arm r, is then given by the relation 


<1^ 


(—a 


(b) 


'I'ho lift-curve slopes and the locations of the aerodynamic 
center vary with the free-sti'cam Afach numlier: hence the 
appropriate values must he used at cacdi flight condition for 
which aeroelastic calculations ai'e made. For airplanes 
designed to oiierate at sul)sonic s])eeds, only the highest 
Mach number attainalile at the highest dynamic [iressure is 
likely to be ci’itical from aeroelastic considerations. For 
airplanes designed to operate at supersonic speeds no such 
general statcnn'iit can be made; however, at a. given altitude 
either the I’cgion of Mach numbers near the ti'ansition from 
the subsonic to the transonic regime or the highest attainable 
Mach nund)ei‘ is lil<ely to be critical as far as the aeroelastic 
|)henomena considered in this report are concerned. (Sec 
fig. 4 of ref. 1 and fig. o of ref. 2, for instance.) 

d'he airspeeds at which the vai'ious aeroelaslic phenomena 
are of interest enter the calculations in tbe form of the cor- 
responding dynamic pressures, d'hese dynamic pri'ssui'cs, 
in turn, are e.\|)resseil in dimensionless foi-m as 


or 


(/ f ( tcA'/' cos A 
1T4 {G^JYr 


(7) 



sin A 

(Ef); 


(8) 


'I'he parameter (/* is useful in the analysis of unswepl wings, 
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for mIucIi loi'sional doformations arc predominant; the 
liai-ameler q is jiriinarily useful for highly swept wings, for 
wliicli liending deformations ai'o predominant. In general, 
(he pammeler q* is used in this report unless ty is zero. The 
ratio of tliese parameters 


TABLE 1.— DEFIXITIO.XS OF THE FACTORS „ TO 

ir„ 


k 


q s, 

q*~f\Cr (Kl), 


tan A 


( 9 ) 


is independent of the dynamic ])ressure and depends only on 
geometric and structural parameters. This ratio is veiy 
us(‘ful for analyzing the aeroelastic behavior of swept wings. 

Structural parameters. — For the purposes of an aeroelastic 
analysis, the wing structure is characterized by the location 
of the elastic axis and the magnitude and distribution of the 
bending and torsional stiffnesses (El and GJ), as well as by 
lh(> magnitude of the rigid-body rotations imparted to the 
wing by its root. 

'I'lie elastic axis is usually defined as the locus of points at 
which normal loads can he applied without causing the wings 
to twist. Such a locus does not generally exist for practical 
wings; however, foi' unswept wings without cutouts an axis 
can he determined which approximately satisfies this con- 
dition. Similarly, for swept wings without cutouts an elastic 
axis can be defined for the outboard part of the wing if the 
wing is considered to be clam])ed along some such line as the 
e(fectiv(> root shown in figui’e 1. In most aeroelastic calcula- 
tions the locus of shear centers for both swept and unswept 
wings is assumed to be the elastic axis. If the structure has 
large cutouts which result in sudden changes in the stiff- 
ne.sses and in the shear center along the span, the charts of 
this rc])oi-t cannot be used except in a cpialitative sense. 

The magnitude and the spanwise distributions of the 
bcjxling and torsional stiffnesses enter aeroelastic calcula- 
tions by means of the charts and approximate formulas in 
diffcM-ent ways. The magnitudes, as characterized by the 
values of the stiffnesses at the effective root, have to bo 
known in oi-dei- to jjerform any calcidations; the distributions 
ai'c imjdicit in the chaids. The root stifl'nesses, if not knomi 
ollu'iavise, can 1)(> estimated either from exjierience with 
similai- designs, from the results of tlie constant-stress concept 
outlined in ajipendix B, or from a combination of the two. 

d'he r(‘quired bending stiffness at the root (EI)r is propor- 
tional to the design load factor, the wing loading, the wing 
thickne.ss ratio, the fourth power of the root chord, the 
sipiare of the swept-span asiiect ratio (A/vos-X). and the 
I'alio of the modidus of elasticity to the allowable bending 
stress and depends on the taper ratio and on the detailed 
design of the wing (see aj)])endix B). By means of this 
I’clation the bending stifl'ness of one wing can be estimated 
fi'om that of a reasonably similar whig. Or. with the con- 
stants of pi'oportionality t]-,o and Fu, given in table 1 and 
aiipendix B, res|ieetively, which take mto account some of 
the detailed design iiarameters as well as the taper ratio, 
the stiffness can be estimated directly. Howevei-. in view of 
the fact that these constants have been derived on the basis 
of a highly idealized structure and loading condition they 
must he used with caution. The ratio of the root benduig 


VI 


1 


Vo = 


56 = 


1 + Margin of safety 

_Ordinalc of most highly stressed element 


One-half of wing thickness 

.Actual skin thickne.ss of most highly stre.sscd element 
Equivalent skin thickne.ss of most highly stressed element 


_ Cross-sectional area of (assumed) single torsion cell 
' Chord X Wing thickness 

_Pcriniclcr of torsion cell 
® Twice the chord 

Effective perimeter p of torsion cell 
* Actual perimeter of torsion cell 

_ Width of equivalent sheet 
“ Chord 


_ -Average ordinate of upper skin 
.Ma.ximum ordinate of upiierskin 

Equivalent thickness of lower skin 
Equivalent thickness of upper skin 

_ Ma.vimum ordinate of lower skin 
One-half of wing t hickness 

_ Average ordinate of lower skin 
Afaximum ordinate of lower skin 

’)i5 = tti(i(>>5-’)n--r ttu’tiritn'-) 

-Allowable torsion st ress 
-Allowable bending stress 

cb' 

’"■^Scos.I 

m.=Vi. Q) 

1 

— - 


mvh 

n:o= ^.r 

V4Vi7‘ni'i' 

_ Equivalent perimeter of torsioji cell 
Actual perimeter of Torsion cell 

stiffness to the root torsional stiffness can be estimated 
b3' means of cejuations given in appendix B or. preferably, 
from experience with structures similar to that under 
consideration. 

The spanwise stiffness distributions need not be kno\^^l in 
detail in order to use the charts and approximate formulas. 
If the wing is solid or nearh' solid or if its cross sections are 
geoinctricallv similar at all points, the charts for stiffness 
distributions proportional to the fourth power of the chord 
arc used. If the wing does not have large cutouts and is 
designed for a constant stress level, the charts for the stiff- 
ness distributions associated with constant stress are used. 


j-i 2 
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The use of tliese eharts (ends to ovei'estimate aeroelastie 
effects to some extent because, although actual wings are 
designed for constant stress over most of the span, the por- 
tions near their tips are designed on the basis of other con- 
siderations, such as handling loads or minimum standard 
sheet thicknesses; thercfoi'c. wings tend to be stift'er near the 
tip than they woidd be if desigjied on the basis of constant 
stress throughout. This difference in stiffness is particu- 
larly large if the taper ratio is zero. 

If the wing contains large cutouts or if, for any other 
reason, the wing stiffness distribution is kno^\^l to be sub- 
stantially different from a constant-stress type, the eharts 
can be used to furnish semiqualitative residts for the various 
aeroelastie phenomena by using fictitious stiffnesses, pro- 
vided the actual stiffness distribution is known at least 
approximately. The root stiffnesses of these fictitious dis- 
tributions may be assumed to lie the ones that give rise to 
twist or bending angles at the tip which are the same as 
those of the actual wing if the bending moments or torques 
vary as the square of the distance from the tip. For con- 
venience, the spanwise distribution of these fictitious stiff- 
nesses may then be assumed to be proportional to the fourth 
power of the clioi'd. On the basis of these assumptions. 



where the subsci'ipt e refers to the fictitious stiffness, and 
where the integral represents the moment of ineidia of the 

area under the function ^ plotted against s* about the 

point s*=l. The fictitious torsional stiffness at the root 
can be obtained in the same mannei'. The aeroelastie 
phenomena can then be estimated by use of these fictitious 
root stiffnesses and the charts for stiffness distributions 
proportional to the fourth power of the chord. 

In the derivation of the charts the wing is considered to 
be clamped at the effective root perpendicular to the elastic 
axis. From the data and anabases pivscnted in references 
1, 2, 4, and 5 a satisfactory location of the effective root 
appears to be at the intersection of the elastic axis and the 
side of the fuselage. 

If the rotations at this effective root are known as a 
result of deflection tests or a detailed anah^sis such as that 
of reference 5, the root twist due to toique and (he roo( 
bending due to bending moment may be taken into account 
by moving the effective root inboard by the distance 

(GJ)r (11) 

or 

( 12 ) 

where v?,’’ is the angle of twist at the root due to a root 
torque 7,, and where is the deflection slope at the 
efl’ective i-oot due to a bending moment ff/,. Since the 


distances A.v'’ and A.v‘' may differ from eacli other, some 
compromise between the two must be made; for unswept 
wings the use of A.s-' appears to be indicated, whereas for 
highly swept wings the use of A.s’’ is more ajtpropriate. 


PRELIMINARY SURVEY OF AEROELASTIC DEHAVIOR 


The information contained in some of the subsequent 
sections has been summarized in figure 2 for the purpose of 
ascertaining in advance of more detailed estimates, if desired, 
whether the aeroelastie phenomena considered heiein are 
likely to affect the design of the wing structure. This pre- 
liminaiy survey is not essential to any of the further calcu- 
lations but may show them to be unnecessaiy in some eases. 

The charts of figures 2 (a) to 2 (d) pertain to wings of 
taper ratios 0, 0.2, 0..5, and 1.0 and constitute plots of the 
(lynamic-i)ressuro parameter q* defined by equation (7) 
against the swec]) parameter k defined by equation (9). 
These parameters contain the root stiff’nesses {GJ)r and 
{EI)^; if, when a preliminary survey of aeroelastie effects is 
to be made, no information whatever concerning the wing 
stiffness is available, the following relations for q* and k 
may be used: 


El Cl (1 -|- X)' cos .S. 


1843.. S Cr 


(13) 


, 1 d- X G 

2^ E 


yi^A 

4ioCi 


rit, tan .V 


(14) 


whore is a root-stiffness parameter defined in appendi.x 
B and plotted in figui'c 3, and where and arc defined by 


and 




VAVsVfiViiiVn 


(15a) 




VbVsVtVti 


(15b) 


in terms of some of the factors defined in table I. 

Figure 2 (e) pertains to wings for which the moment arm 
C( is zero and, hence, k is infinite; witli (he degree of appro.x- 
imation involved in the use of figures 2 (a) to 2 (d), figure 
2 (e) can be used for wings witli |^-] > 25. This figure con- 
sists in a plot of the dynamic-pressure parameter q, defined 
by equation (8), against the taper ratio X. If no informa- 
tion is available concerning the root bending stift'ness (ET)r 
contained in ~q, the following relation may be used: 


(1 +X)^ /1a sin A 


E a IF hr a<i 


q 


3(5864 ^ 

Cb S Cr 7I„ 


•I a T.i 

— Vl9i r 


(16) 


The various lines of the eharts of figure 2 designate the 
conditions at which a wing designed on the basis of strength 
considerations alone is likely to encounter divergence or s])an- 
wise shifts of the center of pressure by various amounts; 
positive shifts are tliosc toward the lip. Tliese spanwise 
shifts furnish an estimate of the increase in root bending 




iwfmm 

^ m/Mt~ 

. 
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-8 -4 0 4 8 12 16 

Sweep poromeler, k 

(a) Wiiifrs of taper ratio 0. 

KiGunK 2.— Charts for a prcliiiiinary survey of aeroelastic plioiioineiia. 


momt'iiL due lo ncfoolastic action and an estimate of the 
shift in n'ino; aerodynamic center, since 

^a = sin A (17) 

e.uvic 

Inasmuch as the parameters 5* and \ contaiti the dynamic 
|)rcssure, the negative values of q* showm in figures 2 (a) to 
2 (d) may require some ex])lanation. The four quadrants of 
(■acli of tlie charts of figures 2 (a) to 2 (d) maj' be characterized 
for practical purposes as follows: 


Qimdr;iiit 

Swfoj) 

1 

Divorgcnco j 

; 1 

Shin in If > 

1 

1 Hack 1 

1 Fosilive 

Impossible beyoiul j 
a certain sweep 

j Inboard beyond a 1 
1 certain sweep i 

1 

1 Forward 

I’osilivc 

Likely 

Outboard 


Hack 

Nojrativc 

Impossible 

Inboard 1 

4 

1 

1 Forwar<l 

1 N'opativc 

' Possible be.vomi a 
certain sweep 

j Outboard beyond :i j 
certain sweep 


For uirswe])! wings k is approximately equal to zero, and flie 
aeroelastic phenomena referred to in the charts of figures 
2 (a) to 2 (d) arc similar to those of swept wings defined bj- 
points in quadrant 2 if 5* is positive and by points in quad- 
1‘ant .3 if 2* is negative. In other words, the aeroelastic 


plienomena of unswept wings are similar to those of swept- 
forward wings if c\ is positive and to those of sweptliack 
wings if C\ is negative. The aerodyjianaic-center shift asso- 
ciated with tlie shift in the lateral center of pressure Tj or in 
the spanwise center of pressure ? is always forward, except 
for small positive values of k (associated with sweep angles 
smaller than a certain value), in both quadrants 1 and 4 . 

The significance of negative values of 2* is that c 1 is nega- 
tive, rather than that 2 is negative. A negative value of ci 
may be obtained at supersonic speeds, but under most con- 
ditions Cl is likely to be positive. Similarly, in figure 2 (e) 
a negative value of 2 implies that .V is negative (that is. that 
the wing is swept forward), whereas a positive value implies 
that the wing is swept back. 

In usmg figures 2 (a) to 2 (d), estimates must be made of 
either the root stift'nesses (in conjunction with eqs. ( 7 ) and 
( 9 )) or of the effectiveness factors and t/o (for use in eqs. 
( 13 ) and ( 14 )). The factor Fr is obtained from figure 3 for 
the largest value of Ci likel.y to be encountered at the design 
load factor and for the given taper ratio X. The parameter 
2* is calculated for the combination of dynamic pressure 2, 
lift-curve slope , and moment arm Ci which is likely to 

be critical from an aeroelastic point of view. For an unswept 
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(h) Wings of lapcr liilio 0.2. 
i'lor UK 2. — ('oiiuniu'd. 


wing: llu' coinhinatiou for whicli the piodiict a 

inaxiinum is lilcaly to l)c critical; for a swept wing: the coml)i- 
nation for which qC\ is a ina.xiniuni is liUelv to !)(> critical. 

The parameter k is then calculated foi' the same value of ci. 

The values of q* and k define a point on one of the charts 
of fig:ures 2 (a) to 2 (d) (whicliever is closest to the actual 
taper ratio). If the shift in spanwise center of pressure (and 
any associated shift in the aei’odynamic center) at tliat point 
is small and, in the case of an unswept or a sweptforward 
wing:, if the ahsoluti' value of the ratio of the value of q* at 
that point to the value of r/* at divergence foi’ the value of 
k at tliat point is small, the static aeroelastic phenomena 
diseussed in this report probably need not he taken into 
aceount in designing the wing structure. On the other hand, 
if tlie point on the chart indicates tin' liktdihood of significant 


aeroelastic ed'ects on the spanwise center of pressure or the 
possibility of an approach to the divergence condition, further 
calculations are desirable. The charts of tliis report may be 
used for the preliminary calculations; once the strueture has 
been designed, inoi-e refined methods such as that of reference 
I may he usi'd. 

If the inomeTit arm e, is so small or the angle of sweep so 
large that the iiaramete]- k e.xceeds the range covered by 
figures 2 (a) to 2 (d), the chart of figure 2 (e) may be used for 
the purpose of a ])reliminaiy aei oelastic appraisal of the gis'cn 
wing. In this figure only the jiarameter g is recpiired, since 
k is considered to he infinite. The parameter g can be ob- 
tained from erpiation (Hi). The analysis then proceeds in the 
same manner as for figures 2 (a) to 2 (d) using, in the parameter 
g, the condition for maximum qC[,^ . 
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(c) IVings of taper ratio 0.5: 
Fioi'nE 2. — Continued. 





10 


REPORT 1 HO — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



(cl) Wings of taper ratio 1.0. 
Figure 2. — Continued. 
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(o) Wings with inoincnl arm (>i=0 (or \l; ^'2.5'. 
Ffoi re 2. — Concluded. 





0 .2 .4 .6 .8 10 


Taper rolio.X 


Figure 3. — The root-stiffness function F^, 
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CALCULATION OF THE VARIOUS AEROELASTIC PHENOMENA 


Dynamic pressure at divergence. — Tlic solutions for the 
divergence speed obtained by tiie direct method in reference 
2 and those obtained by the nnincrical matri.x method given 
in appendix A can be summarized b}- approximate formulas 
which give the dimensionless parameters suid Qd (the 
values of the parameters defined in eqs. (7) and (8) that 
correspond to the value of the dynamic iircssui-e q at diver- 
gence) in terms of their ratio k defined by equation (9). 

These approximate formulas are 


and 


q'D = 


<JD=- 


K, 

-K.k 

(18) 

lu 

r. 1 

(19) 


Wlien the angle of sweep is zero, equation (18) reduces to 
q*D=Ki, and when the moment arm c, is zero, as it may be 

in supersonic How, equation (19) reduces to qo = — The 


constants and K-> are given in table 2 for wings with taper 
ratios of 0, 0.2, 0.5, and 1 .0 for both tyjies of stiffness distribu- 
tions; the parameter for unsweiit wings and the parameter 
qo for swept wings with Ci = 0 are plotted against the taper 
ratio X in figures 4 (a) and 4 (b), respectively. 

With the values of q*o or qo given by equations (18) and 
(19) and the definitions of these two ])arametcrs given by 
equations (7) and (8), the values of q required for divergence 
may be determined. If desired, the corresiionding airspeed 
may be detei-mined from the relation 




/ 

A' p/2 


The value of qu is often negative for sweptbaek wings, and, 
since a negative dynamic pi'ossure does not correspond (o any 
real speed, these wings cannot diverge. These negative 
values of qo, nonetheless, are useful as reference values in 


TABLE 2.— VALUES OF THE COEFFICTEXTS /C, AX'D AX 





A'i 


'Paper 

rano, 

X 

GJ and Ff 

(hlc)t 

(/i./cv 

, lly an.i- 
H\' mall i\ lylic itilc- 
inlC^'ralion pr.it ion 
(I of. 2) 

1 

1 

Hy .ina- 
Hymatn.x lylic iiilo- 
iiitepntiioni pnttioii 
1 (ref. 2) 

1.0 

.5 

2 

o’*" 

Vary as 

1.0 1 
1.0 ' 
1.0 1 
J.O 

2. oS 1 2. 47 

2. H'i 2. 74 

2.92 2.K1 

1.45 [ 2.25 

0. .tsi 
' .480 

.51K) 
.520 

0.390 

.497 

.r*14 

1.0 

.5 

o 

0~ 

Given by constani- 
stross criterion 

1.0 , 
1.0 
1.0 

1.0 1 

.5 

, 79.i 

I.2S7 1 

1.S30 

1.440 

, 92S 

. 252 

I . 357 

' .480 

. 620 
.310 


.5 

Increased beyond val- 
ues required by con- 
stant-stress criterion 

1.0 

1.700 I _ 

i 

.39S 

1 



other aeroelastie phenomena. (If the negative value of <7« 
is very low in absolute magnitude, divergence in a higher 
mode associated with a positive qo is conceivable. However, 



(a) Effccc of t.ipcr ratio on q*u for unswept wings. (?*i>=A', in flii.s 

ease.) 

(I>) Effect of taper ratio on qo for Ci = 0. ^?o = — ^ in this caso.^ 

Ficcre 4. — Etfect of taper ratio on the dyiiainie pro.s.siire at tlivorgeme. 



CHARTS AND APPROXIMATE FORMULAS FOR ESTIMATION OF AEROELASTIC EFFECTS ON WING LOADING 


13 


for Llie orclinaiy wings wMth straight leading and trailing 
edges and substantially straight elastic axes to which the 
analysis of this report is applicable, the magnitude of the 
lowest positive qo is always much larger than the absolute 
magnitude of the negative qo associated with the first mode. 
TTiglier-modc divergence of a practical wing for which a neg- 
ative qo is indicated by the anal^'sis of this report is therefore 
vciy unlikely.) 

The values of the constants Ki and Kn given in reference 2 
difi'er somewhat from the corresponding values resulting from 
ihe matrix solution in appendix A. The matrix results are 
probably more significant because tliej^ are based upon more 
i-ealistie aerodynamic assumptions; the and K> values in 
reference 2 tend to give conservative results. 

The value of q^ calculated for any given value of q*o or 
7 b depends on the value of the effective lift-curve slope Cl^ 

or c, and, hence, on the Mach number. As suggested in 

references 1 and 2, the value of q^ calculated at various 
.Mach numbers may be plotted against Mach number. If 
lines of the actual dynamic pressure at several altitudes as a 
function of Mach number are drawn on the same plot, an 
intersection of the divergence line with one of the lines of 
actual dynamic pressure designates possible divergence at 
that value of dynamic pressure, Mach number, and altitude. 
If this plot is on log-log coordinates, the lines of actual dy- 
namic pressure are straight and the ratio of the dynamic 
jiressure at divergence to the actual dynamic pressure at a 
given Mach number and altitude can be scaled oft' dil•ectlJ^ 
(Sec refs. 1 and 2.) 

Spanwise angle-of-attack distributions. — In appendix A, 
an apiiroximate expression is determined for the change in 
angle of attack duo to wing ftexihility. For the additional- 
tyjic anglc-of-attack distribution (a^ is constant) the angle 
of attack due to structural deformation a, is given by 

(/i+F.A/0 (20) 

K Y i- 

Qd 

The functions /i and A/i of the spanwise coordinate s* and 
the function F, of the parameter k are given in figure 5 for 
swept wings with taper ratios of 0,2, 0.5, and 1.0 and with 
tlie two different types of stiffness distributions. For wings 
with zero taper ratio the structural deformation cannot be 
obtained from ecpiation (20), as is pointed out in appendix 

A. However, the ratio k ~ as a function of the spanwise 

coordinate s* is shown in figure 6 for the two different stiff- 
ness distributions, several values of q/qo, and several values 
of the parameter k. 

The spanwise distribution of due to a linear twist 
which may be either symmetric or antis 3 TOmetric, 

is ajuiroximately 


where the functions fj, A/a, and F> are given in figure 7 for 
wings of taper ratios 0.2, 0.5, and 1.0. The angle-of-attaek 
ratio is shown in figure 8 as a function of s*, qlqo, and k for 
wings of zero taper ratio. 

The results of equations (20) and (21) maA* be superim- 
posed. .For example, if the spanwise distribution of a, due 
to rolling is to be found, these equations must be added in 
such pro])ortion that 




pb 


But 


~b b 


* 


so that 



1-^+^ 

b^b 



pb 

w 


( 22 ) 


where pbj2V is the angle of attack at the tip due to roll. 
The distribution of due to roll is then 


ft= ^ [( (i - t) j (/.+f.i/.)] (23) 


Spanwise lift distributions. — If desii'ed, the lift distribu- 
tions can bo obtained for the angle-of-attack distributions 
given in the preceding section by one of the commonly used 
methods of calculating spanwise lift distributions, such as 
that of reference 6 or by the approximate method of reference 
7. However, the following method is simpler and, within 
the approximation of the present report, as accurate pro- 
vided the rigid-wing (7=0) loading lo is calculated by the 
methods of refeience 6 or 7 or is obtained from the charts of 
reference 8 or, in general, b_v an accurate anah’tical method. 

Within the framework of the assumptions made in the 
anah'sis the lift per inch of span is proportional to the local 
angle of attack, so that 

I a, 

j-==l-|-K~ (24a) 

lo 


for geometrical angles of attack which are constant along 
the span, and 


*0 Ofg 


(24b) 


for geometrical angles of attack due to linear twist, where 
and as/ccg are obtained as indicated in the preceding 

«JT ' 

section. 

If no better approximation is available for the loading U 
(which is likelj- to be the case at transonic speeds), it may 
be estimated for geometric angles of attack which are con- 
stant along the span from the relation 


lo~CLj^qcig 


(25a) 


(/•„-!-/?’., a A) 
1 <L" 

Qd 


( 21 ) 


and for all other geometric angles of attack from the relation 

lo^Cj.^cqag (25b) 


2liS8L>5 .‘>4 - - 



f, and A/, 
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Dimensionless distonce along the span, 5* 
(e) Stiffnesses proportional to c‘. 


Ftci RE 5. — Concliiflotl. 
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(a) Stiffnesses proiiortional to c*. (b) Stiffnesses given by constant-stress criterion. 

Figure S. — Tlie angle-of-attack distributions for linearly varying angles of attack for wings of taper ratio 0. 
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Lift and moment coefRcients. — The wing lift coefficient 
the wing-TOOt bending-moment eoefficient Cb, and the 
wing-root twisting-moment coefficient Ct may be obtained 
in terms of their respective rigid-wing values l)v means of 
the following approximate expressions; 




l-f (l-r) 
1--^ 


Cb 


(In 



1 -^ 


(Id 

Cr ' 

1 ^ (1 — tv) 

(Id 




_JL 

(Jn 


(26) 


(27) 


(28) 


where tlie coefficients r, p, and r depend on the type of 
loading. The sidiseript 1 is used for additional-type angle- 
of-attack distributions and the subscript 2, for linear-twist- 
type angle-of-attack distriliutions. 'f'he coefficients r,, ni, 
and Ti are given in figure 9 as functions of the parameter k 
for wings of tapcj- ratios 0.2, 0..5, and 1.0. For M'ings with 
taper ratio zero the ratios of the lift, bending-moment, and 
twisting-moment coefficients to theii- respective rigid-wing 
values are given in figui'C 10 as a function of qlqn for several 
values of the parametei- k. The values of v-,, 1 ^ 2 . and are 
given in figure 11 for wings of taper ratios t).2, 0.5, and 1.0, 
and ratios of the lift, hending-moment. and twisting-moment 
coefficients are given in figure 12 for wings of zero taper ratio. 
'I'lie additional-twist and linear-twist results of etiuations 
(26) to (28) may he superimposed in the same way as those 
of equations (20) and (21). 

The wing rolling-moment coefficient C, is defined as the 

t- ‘ir 

rolling moment of the loads on both wings about the fuselage 
center line divided by qSb. 'riicrefore. 


C,. 


2.1/,. cos .V-r27Vsin — 


(29) 


'I'he angle-of-attack distril)ution due to rolling given in 
equation (22) must he used in finding tlu' values of .!/„ T^, 
and L,n in equation (29). 

Spanwise centers of pressure and aerodynamic centers. — 
'I'he spanwise loca(ion of the c('nter of ]n-essure is given by 
(he distance 


_ b 



(30) 


( n 


( 31 ) 


(Inasmuch as 17* is equal to .v* by virtue of the definitions of 
those dimensionless (|uantilies (see also fig. 1), eq. (31) can 
he considered to he an expression for rj* rather than -v*, if 
desired.) With the values of the hending-moment and lift 
coefficients given in the preceding section, the ratio of i to 
its value for the rigid wing may he calculated from ('ilher of 
the ec| nations 



All 




_ 

1 ^(1-;-) 


(Id 


(32) 


where the shift in spanwise center of pressure Ail is defined 
as « — So, and where tn and vi are used for constant geometrical 
angles of attack and m and r..>, for linearly varying geometri- 
cal angles of attack. 

'Phe shift due to aeroelastic action of the longitudinal 
position of the center of ])ressure associated with a given 
shift of the spanwise center of jiressure is eciual to sin A As. 
'Phe shift in aerodynamic center (]}ositive when I'earward, or 
stabilizing) can conseciuently he calculated l>y substituting 
into ecpiation (17) the values of As obtained from ecpiation 
(32) with values of g, and v^. 

Inertia effects. — Xo charts are presented in this report for 
th(‘ effects of inei'tia on quasi-static aeroelastic phenomena — 
that is. aeroelastic phenomena associated with flight at 
constant acceleration; the manner in which mass is distril)- 
uted varies so widely among different wings that in-ejiaration 
of a generally ap|)licahle set of charts for inertia effects 
appears to he imin-actical at jiresent. Furthermore, e.xcept 
for flying wings, the wing deformations due to inertia loads 
are small comjiared with those due to aerodynamic loads, 
the two types of loads being in about the same ratio as the 
wing weight to the weight of the entire airplaiu'. If desired, 
liowever. inertia effects and the aeroelastic increment in 
these effects can he estimated in the manner described in 
(he following jiaragraphs. 

P'rom the knowm or estimated mass distribution of the 
wing the inertia load f, per inch of span and the inertia torciue 
/,• per inch of span can be calculated for any given normal, 
pitching, or rolling acceleration. Substitution of these loads 
and tonpies for the terms / and lt\C in equations (A3) oi‘ 
(A36) and ('cpiations (,\2) or (.V.35), respectively, yields the 
values of the accnmulated bending moment and torque due 
to the disti'ihuted inertia loads and torques. In turn, sub- 
stitution of these accumulated bending moments and torqties 
in equations (.A4), (.Ao), and (A6), or in ecj nations (A37) 
and (A38) and the matrix equivalent of equation (A6), 
yields the angle-of-attack distrilnition due to the deforma- 
tions caused by the inertia effects associated with the given 
acceleration. 


or the dimensionless distance 
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(c) Stiffnesses proportional to c*. 
Figure 9. — Concluded. 


Asymptotes 
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(a) Stiffnesses proportional to rL (b) Stiffnesses given b.v constant-stress criterion. 

Figure 10. — The lift- and moment-coefficient ratios for con.stant geometric angles of attack for wings of taper ratio 0. 




(a 
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(b) Stiffnesses related to those given by const ant -stress criterion for wings with taper ratio 0.5. 

Figi re II. — Continued, 


Asymptotes 
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(c) Stiffnesses proportional to r‘. 
Figtre 11 .' — Concluded. 


Asymptotes 
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This angle-of-attack distribution can be considered as a 
geometrical angle-of-attack distribution. For the purpose of 
calculating the increment caused bj^ aeroelastic action, this 
distribution can be approximated by a linear-twist angle- 
of-attack distribution with a value at the wing tip which 
is such that the moment about the effective vdng root of 
(.lie area under the linear-twist distribution equals the mo- 
ment of the area under tlie calculated angle-of-attack dis- 
ti-ibution due to inertia effects. (The moment, rather than 
the area, is suggested as a basis of correlation because the 
angles of attack near the wing tip are more important in 
aeroelastic phenomena than those at the wing root.) The 
justification for this rather arbitraiy approximation to the 
angle-of-attack distribution is as follows: As previously 
mentioned, the wing deformations due to inertia loads are 
likely to be small compared with those due to aerodynamic 
loads; furthermore, the correction to be applied to these 
deformations as a result of aeroelastic action is usually 
small compared with these deformations and, hence, is veiy 
small in comparison with the total wing load, so that the 
correction need not be calculated as accuratel}' as the cor- 
rection for aeroelastic effects to the rigid-wing lift distribution- 

Tbe angle of attack due to structural deformation as asso- 
ciated with the linear-twist distribution can then be obtained 
from equation (21) and figure 7 or, if X=0, from figure 8. 
The lift distribution associated with the total angle-of-attack 
distribution due to the deformations caused by the inertia 
effects, including the increment in this angle-of-attack dis- 
tribution produced bj'^ aeroelastic action, can then be found 
fi'om equation (24b), in which ag and lo pertain to the cal- 
culated angle-of-attack distribution due to the inertia effects 
(not the linear approximation to this distribution). This lift 
distribution can be integrated to obtain the lift, bending 
moment, rolling moment, and aerodjmamic-center position 
due to inertia effects, as modified by aeroelastic action. 

The lift and rolling moment calculated in tliis manner 
may then be combined with the lift and rolling moment for 
steady level or rolling flight calculated b^’' the method out- 
lined in the preceding sections. For instance, if the con- 
tributions of the tail and the fuselage to the airplane lift 
can be neglected, the wing lift can be written as 

where is the total normal force per unit load factor 

due to inertia effects, including aeroelastic effects; it is equal 
to — Wv, plus the lift on both vings due to inertia effects, 
as modified by aeroelastic action, per unit load factor and 
is almost always negative. In the preceding equation 

is a wing lift-curve slope which includes static aeroelastic 
effects and is equal to multiplied by the factor on the 
right side of eejuation (26). Then 

1 1 

" ^l44 , i 7dLX 

\V~W^\dnJs 

~144 ir-iF„ 


where 



^ 1F-1F„U«“A 

is a wing lift-curve slope wliich includes static aeroelastic 
effects, inertia effects, and aeroelastic modification of the 
inertia effects. 

ILLUSTRATIVE EXAMPLE 

The parameters of a swept wing, which differs from the 
wing of the illustrative example of reference 1 only in tlie 
width of the fuselage to which it is attached, are given in 
table 3. The values of As'’ and As^ were calculated from 
tile dimensionless root-rotation constants used in the 
example of reference 1, 0^^=0 and (?r^= — 0.25, by means 
of the relations 

As'’= Qt^Ws 

where Wj, as defined in reference l,is the distance along the 
span between the effective root and the innermost complete 
section of the torsion box perpendicular to the elastic axis. 
In the wing of the illustrative examples of the present 

TABLE 3.— PARAMETERS OF WING USED IN 
ILLUSTRATIVE EXAMPLE 


Structural parameters 

e 

0. UA 

{GJ)t, Ib/sq in. 

..S.WXiO» 

(BDr. Ib/sq in 

..9.56X10’ 

OJl(OJ)r 


EIHEDr - 

...=(c/c,)< 

As*, in 

0 




-3.0 

1 ij in. 

V cos / 

215.9 


Geometrical parameters 

A 

A 

A. deir 

.37. 5 

S, sq in 

7. 498 

5, in 

3S7.4 

w/2, in 

.20.0 

672, in - 

173.7 

Ct, in 

102.8 

Ct, In - 

.54.2 

X 

0. 527 

CiiACt iti 

100.4 


1 

1 Aerodjniamic parameters 


1 

Subsonic 

Supersonic 


(Af<0.65) 

(.v=i.s) 


0.25 

0. 425 

ei - - 

0. m 

0.019 



2. 78 

■1.92 


0. 78 

1.00 




1 

1 

Aeroelastic parameters 



Subsonic 

Supersonic 


(.\r< 0 . 6 ,D 

(.\f=1.5) 

k 

7. 76 

79.0 

A'l- - 

2.82 

2.82 

Ki 

0.474 

0.474 

(I*D 

-1.053 

-0. 0774 

qot ib/sq ft.. 

-0,400 

-2,700 

U 

Fig. 5 (c) 

Fig. 5 (c) 



Fig. 5 (c) 

Fig. 5 (c) 

Fi - 

-0.27 
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0. 662 
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1-0.345 — 

1-0.338 — 


QD 

qu 

Cl^ 

1 « 

1 « 


Qd 

qu ( 


1-0.140 — 

l-O.llS — 


Qi> 


•TO 

1-0.345 -2- 

1-0.338 — 

1 

qu 

Qu 

\ 




32 


REPORT 1140 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


report and of reference 1, lUc is 22.4 inches. These relations 
for As*” and Aa-*' ean be obtained from equations (1 1) and (12) 
of the present report in conjunction with the definitions of 
the root-rotation constants given in equations (15a) and 
(15d) of reference 1; in the notation of the present report 
the definitions are 


Q 


Vt 


Wj(GJ)r 




jZIMl 

Wej{EI)r 


The values of r,, t,, and were obtained from figure 9 (e) 
and substituted into equations (26), (27), and (28). Tlie 
wing lift coefficient, wing rolling-moment coefficient, and 
spanwise centcr-of-pressure ratios, as well as the shift in 
aerodynamic center, were calculated bj’’ use of these ap- 
pro.ximate equations in conjunction with equations (17) 
and (29) and arc shown in figure 12 as functions of the 

dvnamic-pressuro ratio — 

— (Id 

DISCUSSION 


The stiffness is assumed to vaiy as the fourth iiower of the 
chord in the e.xample of the present report. 

The subsonic and supersonic values of the parameter k 
were calculated from equation (9). By means of ap- 
propriate values of the constants K\ and K> taken from table 
3, the values of q*o were calculated from equation (IS) and 
included in table 3. From these values of q*D, the sub- 
sonic and supersonic dynamic pressures at divergence were 
found through the use of equation (7) and are given in table 
3. These values of qo vaiy as the reciprocal of the eft’cctive 
lift-curve slope, the corresponding values of Ci being assumed 
to remain constant. 

In order to find the angle-of-attack distribution for 
additional- tj’^pe loadmgs from equation (20), the values of 
Fi and of the functions /i and A/i were taken from figure 
0 (c). The spanwise change in angle of attack is shown in 
figure 13 for different values of the dynamic-pressure ratio. 




Figure 13. — Effect of aeroelastic action on some aerodi iiainic proper- 
ties of the wing u.'^ed in illustrative e.vanqsle. 


LIMITATIONS OF THE CHARTS AND APPROXIMATE FORMULAS 

The charts and the appro.ximate formulas presented in 
this report arc subject to certain limitations as a result of 
the approximations made in the calculations on which they 
are based. These limitations take the form of restrictions 
on the plan form, on the speed regime, and on the wing 
structure. The results obtainable by the use of the charts 
are likely to be unsatisfactory for wings of veiy low aspect 
ratio or veiy large sweep and relatively unsatisfactoiy for 
wings of zero taper ratio. 

Wings of low aspect ratio are ruled out on three counts: 
(1) the e.xtent to which aerodynamic forces are overestimated 
in replacing the wing by one with an effective root and tip 
is larger for wings of low aspect ratio than for wings of 
high aspect ratio, (2) elementaiy beam thcoiy is unsatisfac- 
toiy for calculating the deformations of wings of veiy low 
aspect ratio (because the effects of end constraint, shear lag, 
shear deformation, and bending-torsion interaction are 
more important when the aspect ratio is low), and (3) the 
assumptions made concerning the lift distribution of the 
wing are more nearly true for wings of high than for those 
of low aspect ratio. 

For wings with very large angles of sweep, also, the use of 
an effective root and tip introduces relatively large errors in 
the aerodynamic forces. Furthermore, the root rotations 
neglected in the calculations (bending rotation due to tor- 
sion and twist due to bending) are likety to be important 
for wings with large angles of sweep. 

The aeroelastic analj'sis of wings with zero taper ratio 
entails certain mathematical difficulties which do not arise in 
the case of wings with nonzero taper ratio. The stiffness of 
such wings is zero at the tip and very low near the tip, so that 
the boundary conditions for (p and F given by equations 
(AlOa) to (AlOc) in appendix A are indeterminate. As a 
result of the relatively large values of the reciprocal of the 
stiffness near the tip, the numerical-integration methods used 
in the matrix calculations are less accurate. These difficulties 
also occur in other methods of solving the aeroelastic etiiia- 
tions, such as energy methods. Furthermore, the structural 
behavior near the wing tip is not represented adequately by 
elementary beam theory. Finally, that the aeroelastic results 
calculated for wmgs of zero taper ratio are not as reliable as 
those for other wings is evidenced also by tbe fact tbat they 
do not lend themselves to systematization by means of 
approximate formulas, as do the aeroelastic results calculated 
for other wings. 

As a result of these considerations delta wings are unsuit- 
able for aeroelastic analysis by means of these charts because 
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tlu'v have low aspoct ratios, largo angles of swoop, and zero 
taper ratio. 

In order to use the eharts two aerodwiamic. parainotors 
mtist be kiiOMTi for any given ease: the effective wing lift- 
curve slope and the section aerochTiainic center. From an 
acModynainic iioint of view the charts of tliis re])ort may be 
us(‘d in almost all eases for which these ciuantitios are known. 
'I'be exceptions stem from the fact tfiat the spanwise distribu- 
tion of the lift on the rigid wing is assumed to be proportional 
to the chord, and the distance from the section aeroihmamic 
(■(■nter to the elastic axis (as a fraction of the chord) is 
assumed to be constant along the span. These assumptions 
are not valid for wings with large angles of sweep and wings of 
low asjiect ratio, as ini])lied previously. They are also invalid 
to a greater or lesser extent for most wings in the ti'ansonic 
region. Consequently, even when the lift-curve slope and 
the section aerodynamic center are known, any results 
calculated for transonic speeds must be used with caution. 

Anotnei- aerodynamic assumption implied in the charts is 
that no concentrated aeroclynamic forces, such as those due 
to a tip tank or nacelle, act on the wing. Relatively small 
nacelles in the inlioard half of the span can probably be 
ignored for the purpose of an aeroelastic analysis at the 
prc’liminary design stage. However, large tip tanks cannot 
usually be ignoi'ed even in a preliminary aeroelastic analysis; 
the aeroelastic jihenomena may in such cases be gi'catly 
umh'restiniated by calculations made with the charts of this 
report. 

'I'he assumption concerning the aiiplicability of elementary 
beam theory to the calculation of wing deformations due to 
a(‘ro(‘lastic action serves to restrict the wings that can be 
analyzed by means of the charts to those of moderate or high 
as])ect ratio, as stated previously. Neglect of cbordwise 
bending (elastic camber) effects in the calculations on which 
th(‘ charts are based serves to impose a lower limit on the 
thickness of the wings for which the charts may be used. 
)Vhether this limit is within the region of practical thicknesses 
is (|uestionable, however, d'he divergence tests of reference 
2, which wei'c iieidormed on flat plates of moderately liigh 
as|)(‘ct ratio and with a thickniess of 2.5 ]iercent, showed no 
obvious cbordwise bending effects, although the relativelv 
small differences between the measured and calculated 
divergence sjieeds may have been due in part to such effects. 

As mentioned iireviously, for wings witli tajier ratios 
between 0 and 0.2 the results of aeroelastic calculations are 
likely to be ladativelv unreliable. For taper ratios greater 
than 0.2, the stiffness of actual wings tends to be gi-eater 
near the tip than that given by the constant-stress criterion; 
consequently, any given aeroelastic effect is likely to be some- 
what less than that calculated on the basis of a constant- 
stress stiffness distribution but much larger than that 
ealculatc-d on the basis of a c* distribution. 

If a given structure contains large cutouts which give 
rise to discontinuities in the stiffness distributions, equation 
(10) can bo used to calculate a fictitious root stiffness to bo 
used in conjunction with charts for c'-ty])c stiffness distri- 
butions, jirovided the magnitudes of the discontinuities are 
known or can be estimated. 

Use of the charts of this report is premised on the assump- 
tion that the elastic axis is at an approximately constant 


fraction of the chord. If the location of the elastic axis 
varies somewhat along the span, the tise of an average value 
tends to give satisfactoi-y results for the aeroelastic phe- 
nomena of swept wings; for imswept wings, however, the 
results obtained on the basis of this apjiroximation have to 
be used with caution. If the elastic axis exliibits abrupt 
shifts along the span as a result of largo cutouts or for other 
reasons, the charts should be used only for moderately or 
highly swept wings. This restriction is mitigated to a 
certain extent by the fact that an abrujit shift in the locus 
of shear centers does not necessarily imply an equally large 
or equally abrupt shift in the elastic axis. 


RELATION BETWEEN STRENGTH AND STIFFNESS AS DESIGN CRITERIA 

The strength of a structure is its ability to withstand 
applied loads without failure; the stiffness of a structure is 
its ability to deform relatively little under the applied loads. 
The two terms arc related (a fact which forms the basis of 
the constant-stress type of stiffness distributions used in 
this report) but are not synonymous. The problem of when 
to design for strength anti when to design for stiffness and 
the related problem of how to design a wing for stiffness 
when required to do so have been recognized for a long time. 
Because of the comjilexity of these problems no generallv 
satisfactory solution exists at present, but the charts ]u-e- 
sented herein shed a certain amount of light on the jiroblem 
insofar as stiffness requirements occasioned by the aero- 
elastic jihenomena considered in this report are concerned. 

The charts of figure 2 indicate the extent to which a wing 
is likely to be affected by aeroelastic ))hcnomena — that is. 
how far it is from divergence and how much its spanwise 
center of |uessure is likely to shift as a result of aeroelastic 
action, provided the icing i-)s designed on the basis of strength 
considerations atone. If the margin against divergence is 
too small, or if the spanwise center of pressure and the 
associated shift, in the aerodynamic center are deemed 
e.xcessive, the wing has to be stiffened beyond the amount 
associated with the required strength. The charts of figure 2 
therefore serve to delimit (he regions in which a wing can 
be ilesigned on the basis of strength considerations alone 
and those in which stiffness considerations pretlominate. 
at. least to the e.xtcnt of satisfying the stiffness rec(uirements 
associateil with the aeroelastic jihenomena considered herein. 

The bending moment of inertia required by considerations 
of strength alone for the root section of a wing is directly 
projiortional to the design load nflF— H’u.), to the spanwise 
coordinate of the center of ])resstire, and to the wing thick- 
ness at the root and is inversely proportional to the allowable 
bending stress /'». .-Alternatively, this bending moment 
of inertia may be considered to be proportional to the design 

, a(ir— H’A , , , . 

wing loading ^ > to the square of the wing area, to 

the wing thickness ratio at the root, and to a function of 

1 - 


the taper ratio! which is 


^whic 


2 \ 

2 if strip theory is assumed to 


(1-fX) 

apply^; it is inversely proportional to Fb and independent 

of the aspect ratio. These relations for the bending 
moment of inertia required by considerations of strength 
alone are implied in the following discussion of structural 
requirements imposed by considerations of stiffness. 
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In general, a wing witli a high value of q* (see eqs. (7) and 
(13)) is most likely to be aft’eeted by aeroelasl icily (see fig. 2) 
and, for a given vnhw of q*, swept wings ai'e imich more 
likely to be affected by aei'oelastieity than unswept ones. 
(vSee fig. 2 and eqs. (9) and (14).) Consequently, the follow- 
ing wings are most likely to lie subject to aeroelastie phe- 
nomena, provided they are designed on the basis of strength 
considerations alone; 

(1 ) Wings operating at a high Hying s]ieed or high dynamic 
pressure 

(2) Swept wings 

(3) Thin wings 

(4) Wings designed for a low wing loading 

(5) Unswept and moderately swept wings with an elastic 
a.xis relatively far back on the chord or likely to fly in a 
condition in which the section aerodi'namic centers are 
relatively far forward on the chord 

(G) Wings operating at a Mach number at which the lift- 
curve slope is relatively high 

For given wing loadings and given wing areas, some 
aeroelastie phenomena of wings designed on the basis of 
strength considerations alone aiqiear to lie substantially 
unaffected by changes in the taper I'atio — for instance, the 
spanwise shift of the center of ])ressure and the dynamic 
pressure required for divergence. (In the case of the dynamic 
pressure required for divergence, the jiarameter q*o (fig. 4), 
the root stiffness, and the root chord deciease with increasing 
taper ratio, and the net eft'eet of taper is small.) On the other 
liand, the change in the lift due to aeroelastie action is more 
sensitive to the taper ratio; it is more significant for wings 
with high taper ratio than for wings with low taper ratio. 

The effect of aspect ratio on aeroelastie phenomena tends 
to be small for unswept wings of a given wing area, because 
these phenomena are determined largely by the magnitude, of 
the parameter q*, which is independent of the aspect ratio for 
a given wing area. For the aeroelastie phenomena of liiglily 
swept wings, however, the parameter q is more significant. 
This parameter is proportional to the swept-span aspect- 
ratio for wings of a given area. Consequently, with a given 
wing area, taper ratio, and design wing loading, the aero- 
elastic effects of swept wings tend to be more jironounced for 
wings with high aspect ratio than for those with low aspect 
ratio. This statement is particularly true for the shift of the 
aerodynamic center, because a given spanwise shift of the 
center of pressure results in a much greater chordwise shift 
in the case of a swept wing of high aspect ratio than in the 
case of a swept wing with low aspect ratio. 

STRUCTURAL WEIGHT ASSOCIATED WITH THE REQUIRED STIFFNESS 

AMien a given wing has been shown to be siiliject to un- 
desnably large aeroelastie effects (by means of the charts of 
this report or by any other method), the problem arises how- 
to distribute the additional recpiired stiffness. If, for 
instance, the dynamic pressure on an unswept wing is within 


10 percent of the dniamic pressure required for divergence 
and a margin of 20 percent is desired, an increase of 10 percent 
in the torsional stiffness along the entire span will produce 
the desired result. The (piestion remains, howT'ver, whether 
structural weight can be saved by increasing the stiffness 
more than 10 percent in some places and less in others. 

Some insight into this problem may be gained, at least 
insofar as the aeroelastie phenomena considered herein are 
concerned, from aeroelastie and weight ealculations that have 
been made for a family of somewhat arbitrarily selected 
stiffness distributions which differ from the distribution 
required by the constant-stress criterion. The ratios of the 
local stiffnesses to those associated with constant stress are 
shown at the to]) of figui'e 14. The structural-weight factor 
F„ is showm for two of these stiffness distributions as a 
function of the taper ratio. The function is proportional 
to the weight lUj of the primary load-carrying structure and 
depends on the manner in which the wing stiffness and thick- 
ness are distributed along the sjian. (See appendi.x B.) 

The results of the aeroelastie calculations for wings w-itli 
taper ratio 0.5, constant wing thickness ratio /i/c along the 
span, and tlu'se two stiffness distributions are included in 
table 2 and figures 5 (b), 7 (b), 9 (b), and 11 (b). The 
designation “excess strength" in these figures refers to the 
stiffness distributions increased over the constant-stress 
requirement, as shown in figure 14, with a value of o — 2.0. 
Tlie results of the aeroelastie calculations for the stiffness 
distributions decreased below the constant-stress require- 
ment to a value of ca—0.5 are the same as those foi‘ the 
constant-stress stiffness distributions for w'ings w-ith linearly 

varying wing thickness ratio and ^^— — 0.5. 

Tlie results of the weight calculations and the aeroelastie 
calculations may be combined in several ways. The dynamic 
pressure at divergence, for instance, can be varied by changing 
the bending and torsional stiffnesses uniformly along the 
span, by leaving the stiffnesses at the root unchanged and 
varying the stiffness distribution in a manner similar to that 
indicated at the top of figure 14, or by a combination of the 
processes. A specified dynamic pressure at divergence can 
therefore be obtained as the result of many combinations of 
root stiffnesses and stiffness distributions. Figure 1.5 (a) 
consists in essence of a plot of the structural weights asso- 
ciated with combinations of this type against the tip stifi'ness 
ratio CO for a specified dynamic pressure at divergence. This 
figure indicates that the least weight is associated with 
values of the tip stiffness ratio greater than 1. Similarly, 
figures 15 (b) and 1.5 (c) consist in essence of plots of the. 
structural weights associated with various combinations of 
root stift’nesses and stiffness distributions required for shifts 
of ±10 percent in the spanwise center of pressure at a 
specified dynamic pressure. Figures 15 (b) and 15 (c) also 
indicate that the structural weight is least for values of the 
stiff ness ratio co greater than 1. 
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Figure 14. — The’structural-weight'function 
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0 .5 1.0 1.5 2.0 

Tip stiffness rolio, <j 

(a) Relative weight re(|uirecl for agivcii dyiiamie piesMiie at rlivergence. 


(1)) Relative weight rccjtiirecl for — = 0.9 at a given dynaniie [)re^^'^lre. 

S.J 

(c) Relative weight required for ^=1.1 at a given dynamic’ pressure. 

I'lcfUK 15. — The effeet of the tip stiffness ratio on the strueiural 
weight recpiired for a given divergence dynandc itre.-isiire or given 
shifts in the spanwi.-.e center of pressure of wings with taper ratio 0.5. 

The significance of (igtife 1 o is that, if a given wing designed 
on th.e Itasis of stfeiigth alone needs to Ire slifi'ened for 
aeroelastie reasons, most of llie stiffening material shotdd lie 
added in (lie outboard regions, provided tlie weight of the 
material other than that of the primarv load-carrying struc- 
ture is unalfeeted by the stiffening process. In fact, on the 
basis of aeroelastie considerations alone, weight might be 
saved in some cases by removing material from the root and 
adding material at the tip; needless to say, howevei'. strength 
rec|tiirements would be violated by this iirocedure. Just 
where the material should be added in the outboard I'egions 
cannot be said on the basis of the calculations made for 
figui-e 15, since these calculations assume any modifications 
to the constant-stress stiffness distributions to be made as 
indicated at the top of figure 14. However, it appears 
unlikely that great weight savings can he had by using modi- 
fications which differ substantially from those of fijrure 14. 

SOME REMARKS CONCERNING THE AEROISOCUNIC WING 

The term “aeroisoclinic” refers to wings which deform 
under an aerodynamic load in such a fashion that the angles 


of attack of all sections relative to the free stream remain 
unchanged. vSuch a wing has the advantage that its aero- 
dynamic loads do not change under aeroelastie action either 
in magnitiule or in distribution; its aerodynamic center, for 
instance, is unchanged, and the wing cannot diverge. The 
achievement of such “section aeroisoclinicism” is very difli- 
cult and can be realized ordy by separate variation of the 
bending and torsion stifl'nesses; even so, the aeroisoclinic 
condition obtains for only one type of aerodynamic loading 
condition at one Mach number. However, an overall type 
of aeroisoclinicism in which bending and torsion action tend 
to cancel for the wing as a whole is relatively easy to achieve. 
This overall type has, for practical purposes, the same 
advantages as section aeroisoclinicism, in that the aero- 
elastic phenomena considered in this report tend to be 
negligibly small for sueh a wing. 

As may be seen from figure 2, at a small positive value of 
the parameter Ic the values of the parameter q* for diver- 
gence as well as those for given shifts in the spanwise center 
of pressure tend to infinity. This particular value of k 
represents aeroisoclinic wings in the overall sense; from 
equations (IS) and (19) it may be seen to be the reciprocal 
of the value of K, given in table 2. Hence, fiom the definition 
of k (eq. (9)), 


•s, {GJ). 

e,Cr {ET)r 


tan A = 


1 


(.33) 


with the implieation that the distributions of the stiffness 
are of either the c‘ or the constant-sti’ess type and that A'o 
pertains to eithei- of these types and to the aiipropriate taper 
ratio. Ecjuation (33) indicates that, for a given plan form 
with assigned values of -v,. c,, and A, the dis])Osable jiaram- 
eters for the achievement of aeroisoclinicism are the elastic- 
a.xis location e. which enters into the parameter f|. and the 
{GtT) 

root-stifi'ness ratio the aerodynamic center is not under 

the control of the designei- to any ajipreciable extent. 

A ilecrease in the toi'sional stiffness can sometimes be 
(‘fl'ected without deci-ease in the bending stiffne.ss or impair- 
ment of the strength characteilstics of the wing, and ov(>i-all 
aeroisoclinicism may be achieved in this manner for swept- 
back wings. Or, if aei’oisoclinic conditions are considered 
at the out.set, a wing can be designed with the elastic-axis 
location relatively far back (in the case of a swept back wing) 
or forward (in front of the aerodynamic center, in the case 
of a swept forward wing) in order to achieve aeroisoclinicism. 
However, the fact that only certain types of aeroelastie 
phenomena are considered in this report must be kept in 
mind. Locating the elastic axis far back or decreasing the 
torsional stiffness, for instance, may lead to flutter difficul- 
ties, the solution of which may recpiire exce.ssive mass 
balancing of the wing as a whole. 


RELATION OK THE CHARTS TO DESIGN I’ROCEDUKE 

The first step in the design of a wing structure, once the 
wing geometry anil the overall aiiplane characteristics 
have been deciiled upon, usually consists of a rough appor- 
tioning of structural material along the span in a manner 
intended to satisfy strength requirements approximately. 
At a later stage in the design procedure the structure is 
checked for aeroelastie effects and modified, if necessary. 
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Tlu> niodiliciilions iiio then checked again, and so on, until 
hotli stin'ness and strengtli requirements are met with what 
is believed to he a near-oiitiimnn structure from weight 
considerations. The charts of this re])ort may he tised to 
facilitate the pi'ocedui-e at several stages. 

At the very outset, the preliminaiy-stirvey charts can he 
u.s(‘d to ('stal)lish some ovei'all aeroelastic characteristics of 
the wing structtire that wotild he obtained by designing the 
wing for sti-ength alone. If these characteristics are satis- 
factory, the design of the wing structure can proceed on the 
basis of strength rec|uii-ements alone. The final design can 
then he checked for the aeroelastic effects considered in this 
i-(‘poi't hy means of the eharts eontained herein, and for other 
aeroelastic efl'ects, such as flutter and loss of lateral control, 
hy equally approximate methods. However, if the pre- 
liminary survey indicates that a wing designed oti the basis 
of strength alone wotdtl he tmsatisfactory from consideration 
of aeroelaslicity, sufficient additional stiffness may be incor- 
])orated in the preliminaiy design stage, provided the taper 
ratio does not differ greatly from 0.5 and the wing thickness 
I'atio isconstant along thespan. Forinstance, the preliminary- 
survey charts may indicate a shift in the spanwise 
(■(•nter of jtressure which gives rise to a shift of 4 percent in 
the aerodynamic center; whereas the desired maximum shift 
is ‘2 jtercent, so that the spanwise shift must he reduced to 
50 percent of that indicated on the preliminary-stirvej’^ chart, 
d'he shifts in the spanwise center of pressure for a wing with 
inert'ased stiffness at the tip (the “excess strength” case, for 
which o) = 2.0) and for a wing with decreased stiffness at the 
tip (llie wing with to = 0.5, for which the results of the case of 
(li/c), 




= 0.5 may he used) can then be obtained from figure 


9 (h) and e(|uation (.‘12), in conjunction with the value of the 
dynamic ])resstire at divergence estimated from equation (18) 
or (19). The fact that the wings with co=2.0 and o>=0.5 
liav(- dillVrent dynamic ])ressures at divergence than does 
the constant-stress wing must he kept in mind. 

From the .shifts of the spanwise center of presstire cal- 
culated in this manner the value of w for the desired spanwise 
shift can he obtained hy interijolation and, hence, the approx- 
imate magnification factors to be apjtlied to the stifl’ness 
dislrihution for constant stress can he obtained from the 
cliart at the top of figure 14. Estimates for the other 
aeroelastic characteristics considered in this leport can then 
he obtained for the wing with this modified stiffness distribu- 
tion hy interpolating between the restilts given for these 
characteristics for wings with w=0.5, 1.0. and 2.0; that is, 
for the cases referred to, res])ectively, as 


(/(/c), 

(/i/c), 

(/i'C), 

(hU')r 


= 0.5 


= 1.0 


and 


(hin, 

(hic). 


= 1.0 (excess strength) 


in table 2 and figures 5 (b), 7 (b), 9 (b), and 11 (b). Once 
the stnicttire of such a wing has been designed, the various 
aeroelastic effects considered herein should be checked by a 
more accurate method, such as that of reference 1, and the 
loss of lateral control and the flutter characteristics shotild 
be calctilated. 


CONCLUDING REMARKS 

Charts have been jiresented for the estimation of aero- 
elastic effects on the spanwise lift distribution, lift-curve 
slope, aerodynamic center, and damping in roll of sweiH 
and unswept wings at stibsonic and supersonic speeds. 
Two tyjtes of stiffness distributions have been considered; 
one which consists of a variation of the stiffness with the 
fourth power of the chord and is appropriate for solid wings, 
and one which is based on an idealized constant-stress struc- 
ture and is believed to be more nearly representative of 
actual structures. 

The limitations of these charts are that they do not apply 
to wings with very low aspect ratio or very large angles of 
sweep nor to wings with huge sources of concentrated aero- 
dynamic forces. The cliarts are likely to be less reliable 
for wings with zero taper i-atio than for wings with other 
taper ratios and less reliable when the component of the 
Mach number perpendicular to the leading edge is transonic 
than when this component is either subsonic or supersonic. 
Wings with large discontinuities in the spanwise distribution 
of the bending or toi’sional stiffnesses cannot be analyzed 
directly by use of the charts, but a means of making aj)])roxi- 
mate calculations for such wings has been presented. Xo 
charts have been presented foi’ inertia efl'ects hut a method 
of estimating these effects has been outlined. 

In addition to facilitating the calculation of various static 
aeroelastic phenomena, the charts serve to simplify design 
procedure in many instances, because they can be used 
at the jireliminary design stage to estimate the amount 
of additional material required to stiffen a wing which is 
strong enougli and because they indicate that the best way 
of distributing this additional material is to locate most 
of it near the wing ti]). 

Mso. the charts facilitate the achievement of aeroiso- 
clinic conditions, inasmuch as they serve to define a sim|)le 
relation between the elastic-axis location and the wing 
stifl'ncss ratio which is required to obtain this condition 
for a given plan form. Finally, the cliarts indicate that a 
wing which is strong enough is most likely to be affected by 
aeroelastic phenomena if it is to operate at liigh dynamic 
pressures, if it is thin, if it has a large angle of sweep, if it is 
designed for a low wing loading, if it has an elasti(‘-axis 
location relatively far back on the chord, and if it is to 
ojierate at transonic or high suiiersonic Mach numbers. 


Laxglkv Aeronautical Laroratory, 

X’ational Advisory Co.mmittee for Aeronautics. 
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APPENDIX A 

METHODS OF CALCULATIONS ON WHICH THE CHARTS ARE BASED 


THE AEROELASTIC EQUATIONS 

The methods of calculating aeroelastic phenomena usetl in 
preparing the charts of this report are based on the following 
assumptions: 

(1) Aerod^’namic induction is taken into account l)y apjjly- 
ing an overall correction to strip theory and, when matrix 
integrations are used, by rounding off the resulting load 
distribution at the tip. 

(2) AerodATiamic and elastic forces are based upon the 
assumption of small deflections. 

(3) The wing is clamped at the root perpendicular to a 
straight elastic axis (see fig. 1), and all deformations are con- 
sidered to be given by the elementary theories of bending 
and torsion about an elastic axis. 

In keeping with assumptions (1) and (2), the force per unit 
width on a wing section perpendicular to the elastic axis is 

(Al) 

where a, and ag are, respectively, the angle of attack due to 
structural deformations and the rigid-wing angle of attack, 
in planes parallel to the plane of symmetry. (The. geometri- 
cal angle of attack is considered to be constant along the span 
in equation (Al); in the case of linear twist the coefficient 
is used instead of The torque of this force about 

the elastic axis is IciC for uncambered sections. 

The integral equations for the accumulated torque and 
the bending moment are 

r=J^*'/e,Cf/.s- (A2) 

d/=J (A3) 

and, insofar as assumption (3) holds, the angles of structural 
twist and bending referred to axes parallel or perpendicular 
to the elastic axis are 




'(h 

(A4) 


(Ao) 


The angle of attack due to structural deformations is related 
to (/> and r b3’' the equation 


a^ = <p cos A — r sin A 


(A6) 
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Combinuig equations (Al) to (A(5) gives two simultaneous 
differential equations: 

^ J (i?/ cos A- r sin A)1 (AS) 

These equations are subject to the following boundarv 
conditions: 

Zero twist and bending at the root, 

¥>(0)=0 (Afla) 

r(0) = 0 (A9b) 

Zero toi'que. moment, and shear at the tip, 

=0 (AlOa) 

\ (Ifi = 

=0 (AlOb) 

\ aN /s=Si 

In the following sections, equations (A7) and (AS) arc 
solved explicith' for an untapered wing with constant stiff- 
ness along its span and bv matrix integration for a whig with 
an\' arbitraiy stiffness and chord variation. 

SOLUTIONS FOR UNIFORM WINGS 

Arbitrary geometric angle of attack. — -If the torsional 

stiffness, the bending stiffness, and the chord of the wing 

have constant values {GJ)r. (EI)r. and iv, respectiveh", along 
the. wing span, equations (A7) and (AS) become 

,p" cos A= - (j* ag+{^ cos A-r sin A)J (Al 1) 


r"' sin X— — q |^~ a 

where the differentiation denoted b_v the primes is with 
s 

respect to ^=1—— and the dimensionless parameters q* and 
S/ 


,,-l-(.,5 cos A — r sin A) 


(A12) 


q are defined by 


Cl^ qe,cri<,- cos \ 

lATCGJfr 


EL^qCrh^ sin A 
l44 (EI)r 


(A 13) 
(A 14) 



CHARTS AND APPROXIMATE FORMULAS FOR ESTIMATION OF AEROELASTIC EFFECTS ON WING LOADING 


39 


DidVrc'iUiating equation (All) once with respect to ? and 
combining it with equation (Al2) yields tlie single difl'erential 
equation 


///I , ♦ / ~ I T, S 

(Xs — (/ 


('I'lie factor - is used with a, for tlie sake of eonsistenev, 

K 

des|)itc the fact that a geometrical angle of attack which is 
constant over the span does not have a spanwise derivative.) 
Equation (A15) is subject to the following boundarv 
eonditions; 

From eqtiations (.A9a) and (A91)) 


0!s(l ) = 0 

From equations (AlOa) and (AlOb) 

«/(0) = 0 

From equations (AlOc) and (.All) 


(A16) 


(A 17) 


(0)=-q* + (0) 


wliere functional notation is used, so that, for instance, a, (1) 
means tbe valtic of a, at 1 . 

'I’iie solution of equation (Al5) can be effected veiy readily 
by means of Laplace transforms. The complete solution of 
tins equation is 

a. (?) = ™./3 ($)-//(?) (A19) 

where' tlie integral H(^) is defined as 

II (?) = [?*/. (?- ?.) - 7.L (?- ?.)] (A20) 

The functions/a, /^, and/s are defined by 

.h (?) ^ Cxe -=«+c« (r, cos T?+ Y si'i T?) (A21a) 

/A?) = (<e--«+c«(r5COS7|4-^®sin 7 ?) (A21b) 

m) ^ + (f's cos 7 ?+^^ sin 7 ?) (A2lc) 

wbere —2/3 and fi±iy are the roots of /■^-{-q*r—q = 0 and 

‘ 9^0 + 7 ^' 

C,= l-C, 


3|3"-g7- 


' 9^- + y- 


c,=-c. 


^ _ 3g- + 7- 
' 9|3^-f7^ 


C-= — 

' 9)3-’-i-7- 

c,= -r. 


3 <3 

'9/3--f-7' 


The condition for divergence is that ccj bo finite when otg is 
zero along the entire span. .4s can be seen from eqtiations 
(A 1 9) and (.420), divergence can occur only when 


,L(1)=0 


(.422) 


Thus, for a partietdar value of the parameter 

of q* (or q) at divergence is the one which satisfies eqtiation 
(A22). 

Constant geometric angle of attack. — For the additional- 
typo angle of attack, ag(?) = Constant; 


//(?) = li-.L(?)l- 

K 


n(xA^) + 

<x. .fall) 


(A 2 3) 


(A24) 


For lift distributions based on assumption (1) given at the 
beginning of this appendix, the lift per unit width of span 
may then be written as 

fo .fsf 1 ) 

The wing lift coefficient, the wing-root bending-moment 
coefficient, and the wing-root twisting-moment coefficient are 
given in general by 




X Jo' i 
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_4M. 
^ '^~ gSb 




XlCr 

S /2 6/2 





(A26) 


Cr — 


2Tr 

qScr 




a 


■^Cj. 

s/2 


e\ 





Thou, for the imiforiu M’ing, 



Ka M) + 


/^(l) 

/Al) 


and 





- 9 /^( 1 ) 
“ Mi) 




(A 2 7) 


(A28) 


(A 2 9) 


Linearly varying geometric angle of attack. — For the 
linear-twist-type angle of attack, ag(?) = (1 — the factor 
K is 1 . and 

H(^) = ag{^)—\M0—Mb)\ (A30) 

so that 

h<, a's, 

=|J|J/3(0-/A£) (A31) 

The ratios of the wing lift, wing-root twisting moment, and 
M'ing-root bending-moment coefficients to their rigid-wing 
values are then, on the basis of assumption (1). 


r. 

r 


Cb. 


and 


^9 P cvdS)-ro;?(0 

L ^ *■* 

’ = .3 r [ ^ j 

r„ Jo Jo 


('t ^ i 


as in the preceding section. 


M C\ 

0 iro 


(A 3 2) 


(A 3 3) 


SOI.UTIONS FOR NONUNIFORM WINGS 

Equation (Al) may be written in the matrix notation of 
reference 1 as 


f'vM {«.+?! 


and ecpiations (A2) and (A3) as 
{T}=^i,[7'],/c,c- 


= iT4 

{-l/}=srl//']{/} 
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(A34) 






where the matrix [/'] ]Kuforms an integration of the running 
torque /cic from the tip inboard, and the matrix [//'] per- 
forms a double integration of the running load from the tip 
inboard. These matrices arc derived and given in reference 
1. They are based upon Simpson’s rule with a modification 
at the tip, where the load distribution is assumed to go to 
zero with an infinite slope at the ti]). 

Equations (A4) and (Ao), written in matrix notation, arc 



where the matrix (/]” sei-ves to integrate the accumulated 
toripie or bending moment outboard from the wing root. 
This integrating matrix is based iqion Sinqison's rule without 
the tij) modification and is given in reference 1. 

The substitution of ecpiations (A35), (A3tij, (A37), and 
(A3<S) in th(> matrix ecpiivalent of ecpiation (A6) yields 







where the aeroelastic matrix [A] is defined by 


\A]Miy 




(A40) 


The jiarametcrs q* and q are defined by equations (iM3) and 
(Al4), resjiecl ivc'l.v. and 



IVhen ttg is zero along the entire span, ecpiation (A30) 
becomes 

{ ttj } = cy* [*1 ] { 0 's } ( A4 1 ) 

Consc'quently, for a particular value of k, the value of q* 
at divergence can be found by the iteration of the aeroelastic 
matrix (.1). 
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Kc(iia(ioii (A:59) may be rearranged as follows; 

[lll-r/[A]]{.a, + a,} = {a,} (A42) 

Tli(‘ set of linear siimiltaneoiis equations represented by 
equation (A42) may then be solved for the total angle of 
attack + in terms of the values of a* along the span. 

The integrations in equations (A2o), (A26), and (A27) 
may be ])(>rformed with the first rows of the [/'] and [//'] 
matrices. Thus 


a. 




(K0.5+ax} 


(A43) 


(A44) 




l/'J. 


)1 

{ KUs + } 

' c ^ 

f 

(«*} 


(A45) 


4'he, aei'oelastic characteristics of unifonn wings ivere 
calculated by both the direct method of the preceding section 
and the matrix method given in this section. The values of 
the divergence parameter q*D, calculated b\- the direct 
method, were found to be about 5 percent greater than the 
corresiionding values calculated by the matri.x method. 
This discrepancy can be shown to be almost entirely due to 
the lounding oft’ of the loading of the wing tip in the matrix 
method. The dift’ei'ences between corresponding values of 


CVjC < L 


(\ , Or , 

— ' are negligible. 


,, > and „ 


COMBINATION OF RESULTS 

'I'he foi-ins of the approximate formulas used in combining 
the results of the many computations indicated in the analysis 
were obtained by considei’ing a higldy idealized semirigid 
wing; that is, a wing which is rigid along its entire span but 
can bend and twist at the wing root subject to the restraint of 
a bending and a torsion spring. 

If it is assumed that the two spring constants correspond 

(o and the value of q* at divergence is given by 

Sj Hi 

the siinjile formula 


\-iuk 

where the factors Ki and depend on the taper ratio and 
the spanwisc variation of the stiffness. As shown in reference 
2, this formula serves as a good approximation to the calcu- 
lated values of q*o- 


For the semirigid wing, the ratio of a, to ccg is found to be 
proportional to • Iti order to adapt this expression 

to the flexible wings considered in the present analysis, the 
following approximate expression was found to provide 
satisfactorv correlation: 


— {f+F^n (A47) 

] (j_ 

where / and A/ are functions of the spanwise coordinate x* 
and the wing-chord and stiffness variations; F is a function 
of the parameter k and the wing-chord and stiffness varia- 
tions. The functions/. A/, and F also depend on the tyiie 
of spanwise variation of the geometrical angle of attack, the 
subscripts 1 and 2 being used to distinguish between the two 
types of interest. The accuracy of equation (A47) is 
illustrated in figures 16 and 17. 

If equation (A25) is used for the wing lift coefficient (with 
f replaced by s*) and equation (A47) for the angle-of-attack 
distribution, the wing lift coefficient may be expressed as 


or 


C^=, 


JL 

rjD 


1 — 







1-f (1-r) 
'Id 



(A48) 


(A49) 



FiiiCBB 16. — Comparison of angle-of-attack ratios calculated by the 
matri.x method of appendix .A with those calculated from cc|uation 
(20) for constant geometric angles of attack at various dynamic- 
pressure ratios. X = 0.5; k = 8; stiffnesse.s proportional to c*. 
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0 .2 .4 .6 .8 1.0 

Dimensionless dislonce olong span, 5 * 

Ficvre 17. — Comparison of angle-of-attack ratios calculated by tlie 
matrix metliod of appendix .\ with tliosc calculated from equation 
(20) for constant geometric angle.? of attack and for various values 

of the parameter k. X = 0.o; — = 7 ; stiffnesses proportional to <4. 

9o 4 


where the parameter 


^ */ U ^*7 

Jo Cf ocg^ 


so that V and 7 = 7 -“ are fimctions of k and of the wing-chord and 

stiffness variations and depend on the type of geometrical 
angle-of-attack distribution as well. 

As indicated by equation (A47), within the approximation 
inlierent in that equation, the shape of the spanwisc distribu- 
tion of «j does not vary with dynamic pressure. Therefore, 
to a good approximation, the lateral center of pressure of the 
lift due to (as well as that due to ag) does not change its 
position along the elastic axis when the dynamic pressure 
changes. The following approximate formula for the wing- 
root bending-moment coefficient may therefore be deduced 
from equation (A4S); 


= (A.51) 

j q '"o ''0 

Qd 

where and Sg* are the dimensionless moment arms about 
the effective wing root of the. lifts due to 07 . and ag and are 
defined bx^ 

M, 

— J, * L. 
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Tlie values of v, y., and t are given for the two types of 
gooinelrical angle-of-attack distributions in figures 9 and 11. 

J’igure 18 sliows the apjiroxiinate formulas (A49), (Ao2), 
and (A55) to be in good agreement with more acctirately 
eomputed v^ahies. 

The foregoing apiiro.ximate formulas for the structural 
angle of attack and for the lift, bending-moment, and 
Iwisling-moment coefRcients are not applicable to wings with 
y.om taper ratio. An attempt was made to combine and to 
systematize tlie results calculated for such wings in the 
manner employed for wings with other taper ratios, but the 
aiiiH'O.ximate formulas obtained in this way were found to 
j’ield unreliable results. Consequently they are not pre- 
scMited in this report; instead, the results calculated for the 
wings with zero taper ratio are presented directlj’’ in figures 
G, 8, 10, and 12. 



Fioure 18. — Comparison of lift- and moment-coefficient ratios for 
constant geometric angles of attack calculated by the inatri.': method 
of appendix A with those calculated fiom equations (26), (27), and 
(28). X=0.5; k = 8; stiffnesses proportional to c*. 





APPENDIX B 

STIFFNESS DISTRIBUTION OF CONSTANT-STRESS WINGS 


outline of constant-stress concept 


ASSUMED APPLIED LOADS 


In order lo calculate aeroelastic effects, tlie bending and 
torsional stiffnesses of the wing structure El and <?t/liave to 
be known. These stiffnesses enter the calcuhiLions in two 
ways. The root stiffnesses, as indices of the overall bending 
and torsional stiffnesses, constitute, jjriinarv parameters 
which are required in the use of the charts of this report but 
were not required in the j)re|)aration of the charts. On the 
other hand, the stiffness distributions — that is, the ratios of 
the local stiffnesses along the span to tlie root stiffnesses — are 
secondary parameters which are not rc'quired in detail in the 
use of the charts Imt did have lo be assumed in order to pre- 
pare them. 

In calculations jireliminary to tlie actual design of the 
structure, the bending and torsional stiffnesses of the struc- 
ture are not known; they must be estimated on the basis of 
either past experience or considerations of an idealized struc- 
ture. For the purpose of estimating stiff'ness distributions, 
jiast e-vpcriencc with similar structures is likely to be a. useful 
guide in any specific case hut does not lend itself to generali- 
zation and hence to the preiiaration of generally applicable 
charts. The stiff'ness distributions (other than those which 
vary as the fourth jiower of the chord) used to |jrepare the 
charts of this report have been obtained from considerations 
of an idealized structure, as outlined in this ajqjcndix. 

Basically, the method of this appendix consists in an eff'ort 
to relate the stiffness of a wing to its strength and to estimate 
that strength on the basis of certain assumptions. The 
fundamental assumjitions are that the bending and torsional 
stresses are constant along the span and that the ajiplied 
loading is pro port ional to the local chord. The other assump- 
tions concern the. bending and torsional stresses (aiu.sed by 
this load and their relation to their allowable values. In 
estimating these stresses the structure is assumed to be 
essentially of tlu' tbin-skin, stringer-reinforced shell type. 
Certain eff'ectiveness factors are used — for instance, the ratio 
of the allowable torsional stress to the allowable liending 
stress, or the ratio of the cross-sectional area of the effective 
torsion cell to the product of the chord and the wing thick- 
ness. The root stiff'nesses estimated by the method of this 
aijpendix depend directly on the values of these ratios. The 
stiffness distributions, on the other hand, are largely inde- 
jiendent of these ratios but imjily the assumjition that the 
ratios are approximately constant along the span. Gonse- 
((uently, the constant-stress concejit use.l in this ajrjiendix is 
more likely to furnish useful results for stiffness distrilnitions 
than for the root stiff'nesses, and, because of the type of 
structure assumed, the concept is not ajijilicable to very thin 
M'ings. 
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If the applied normal load is distributed in a manner 
proportional to the chord, that is. if 

!=Kc (Bl) 

the bending moment at any point on the span can be obtained 
by integrating the chord distribution asfollows: 





cdx*(h 




where s'* is the dimensionless distance along the reference 
axis measured from the effective root. Similarly, the total 
normal load on one wing is given by 


cos .V Jo 


If the wing is linearly tapered, so that 
C = Cr [1 — (1 — X)a-*] 
where the taper ratio X is defined by 


(B2) 


X = 


l‘r 


then the ratio of the bending moment at any point of the 
span lo the product of the total normal load and the wing 
semispan less one-half of the fuselage width can be exiire.ssed 
as follows: 

- M (B3) 

cos .\ 


where the functioii /’o of .v* and X is delinea by 


I / 1-1-2X 
3 V T-f ^ 



(B4) 


and shown in figure 19. 

.Similarly, if the moment arm of the normal load ajiiilied to 
the wing at any station is also jiroiiortional lo the chord, the 
constant of pro])ortionality being Ci, the distributed lonpie 
at any station is then c,r/. and the accumulated lorcpie is 


T 


= c,/v 


b'!2 

cos 



which may. in turn, be expressed as 


T 

C\CP 




(B.5) 


(B6) 
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wlii'fp till* funclion /t of .>.•* and X is defined hy 

4 ri +X + X- (2 + X) (1 -X) * , (1 *r\ Q __ /g. 

■^’" 3 L (1+x)' (1+x)- +(id-xy- ' 

and ihe average elioi’d c is defined by 


and the fraction of tbe ^ving weight (including the amount of 
fuel, external stores, and so on used in the critical design 
condition) to the total design gi-oss weight is designated hy 


then equation (B9) may he written as 


The fund ion /7 is shown in figui'e 20. 

'rhe total normal load on one wing P can he estimated 
from the design gi-oss weight and the design load factor of the 
aii‘])lane in the following manner: 

P=~{L^—n W^) (B9) 

If the fraction of the wing lift to the total lift carried hy the 
airplane (including that of the fuselage and tail) is designated 
hv 7)1, so that 

L. 



I .2 .3 .4 .5 .6 .7 .8 

Dimensionless sponwise distance, s* 

I'lccHi; 19 .— Tlie IjoiKling-inoinent function /«. 



where 


43=41 — 42 


With the value of P given hy equation (BlO), equations (B3) 
and (B6) serve to express the local bending and torsional 
moments in terms of known design parameters. 

EFFECTIVE SKIN THICKNESS REQUIRED TO RESIST APPLIED LOADS 

The wing structure has to resist both the applied bending 
moments and the applied torques; in other words, the load- 
carrying members must resist combined axial and shear 
stresses. A relation commonly used in the design of wing 
structures loaded by compressive and shear stresses due to 
bending and torsion moments is 




Dimensionless sponwise distance, 5* 
Fiocre 20 . — Tlic twisting-moment function 


where is the applied bending stress, J, the applied shear 
stress, F„ the allowable (compressive) bending stress, and Ft 
the allowable shear stress. However, a similar relation, 

4+4=1 (Bll) 

I B i T 

is somewliat conservative and much more convenient for the 
present purpose and, consequently, is used as the basis of the 
following development. If the margin of safety is not zero, 
equation (Bll) can be rewritten as 

4+4=41 (B12) 

r B ^ T 

where r^^ is an effectiveness factor which can be expressed in 
terms of the margin of safety (M. S.) as 


’'"~ld-M. S. 

The applied bending stress is 


where 2 is the maximum ordinate on the compression side 
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measured from and normal to the chorthvise principal axis. 
Similarly, the applied shear stress is 

(B,« 

ivhere A is the cross-sectional area of the (assumed) single 
torsion cell and ( is the skin thickness on the compression 
side. Substitution of equations (Bl4) and (Bio) into 
equation (B12) yields 


Iz/ . T F, I \__ 
FrOAlz) "" 


In order to relate the bending and torsion stiffness of the 
wing to the skin thiclcness t or to an equivalent thickness Ic 
which includes the material in the stringers and spar flanges, 
the bending stress is assumed to be carried by a box covered 
with sheet of an effective thickness t^, the webs are assumed 
to carry no bending stress, and the torsion stress is assumed 
to be resisted bj' an equivalent single cell, the two webs of 
which contain all the material of the actual webs. The tor- 
sion and bending moments of inertia mav then be written as 


= 2 cfrL 


(B17) 


I=VwClf ^7)5- J 

= ’7ioC^|j U (BlS) 

where the effectiveness factors 7/5 to iju ni'c defined in table 
1. In the factor 7)9, the effective perimeter J> of the torsion 
cell is the sum of the lengths of skin around the perimeter, 
each weighted bj' the ratio of the thickness of the critically 
stressed element to the thickness of the given length of skin. 

WTen the value of / given by equation (BlS) is substi- 
tuted into equation (BI6), equation (BI6) may be written as 




(h 

^^2 

2Azrj, 


Bj' making use of equations (B3), (B5), and (BlO), as well 
as the effectiveness factors 7)19 to riig defined in table 1, 
equation (Bl9) can be written as 


» IB b '/2 J\ 
ViVis Fg cos -t h/ 




The factor/s is defined in terms of the factors/e and/7 given 
b3’' equations (B4) and (B7) as 


r -ih 

.' S O / 


I I V, 

1 . 0 - 

.5 

— : 0 ,- 


I I I I I 
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I'lru RK 21. — The nioiiicat-ratio function /j. 


BENDING AND TORSIONAL STIFFNESSES 

Sulistitution of the value of U given by equation (B20) 
into equation (BlS) jdelds an expression for the bending 
moment of inertia I or for the liending stiffness El at any 
point along the span. The value of this stiffness at any 
point on tlie span inaA'^ be divided ii.v the stiffness at tlie 
wing root (EI)^. This ratio can then be expressed as 

j - I 7; 1 9C 1 

El /e /s Ax 


where 


{El), 

Ba 


,• I 4 lilt' 1 

f,==- A 

' ** ‘•k f'r r 


The function fg is plotted in figure 22. The value {El ), mav 
be, obtained from equations (BlS) and (B20) as 

iEl),= T "J r ,= ( t) F, (x, (B 24 ) 


wiierc 


ir fi^r f , I 


The function /g is shown in figure 21. 


The function F, is shown in figure 3 as a function of X with 
as the parameter. 

B A 

Similarly the torsional stifl'ness GJ ma.^^ be obtained b.v 
substituting the value of I, given bj' equation (B20) into 
ecpiation (Bl7). However, from ecpiations (Bl7) and (ElS) 
the ratio of the torsional stiffness to the bending stiffness 
max' be obtained in the form 

(B 26 ) 
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This equation shows that the ratio GJjEI is constant along 
(lie span witliin the framework of the constant-stress con- 
cept. Equation (E26) ina}q therefore, be interpreted as an 
expression for the value GJjEI at the wing root, that is, for 
the value {GJ)rl(EI),. The torsional stiffness at any other 
jioint on llie span can then be obtained from equation (B 22 ), 


since 


GjJ _ El 
{GJ), {El), 


(B27) 


liecause GJjEI is constant over the span. 

The stiffness ratios GJj{GJ), and EIj{EI), can be ob- 
tained directl}'^ from figure 22 when the thickness ratio hjc 
of the wing is constant along the span; if the thickness ratio 
is not constant, the factor /g obtained from figure 22 must 

lb jC 

I)c multiplied by the ratio station to obtain the 

stiffness ratio at that station. As may be seen from figure 
22 , the function /g does not vary much with the parameter 


this ])arameter represents the additional amount of 

^1a 

skin thickness required to cany the torque (see eqs. (Bl9) 
and (B20)), and this additional thickness is small for most 
conventional wing structures. Consequent!}'’, an average 


value of —',— = 0.03 was used to obtain the stiffness distri- 

A A 

butions used in the acroelastic calculations on which the 
charts of this report are based. 

Equation (B22) shows that, once a value has been assumed 

for the term^^^^j the stiffness ratios EIj{EI), and GJj{GJ), 

are independent of the effectiveness factors used in this 
analysis. Therefore, specific values of these parameters need 
not be known in order to estimate the stiffness distributions, 
hut one of the assumptions on which equation (B 22 ) is based 


is that whatever values the effectiveness factors have are 
nearly constant along the span. In order to estimate the 
value of {EI)„ however, these factors must be known, since 
they enter directly into equation (B24). Estimates of [El), 
and {GJ)r obtained in this manner are, therefore, subject to 
all the limitations imposed by the approximations of the 
constant-stress concept. Hence, some judgment must be 
exercised in using these estimates, and, if possible, they should 
be modified in the light of expei-ience. 

STRUCTURAL WEIGHT ASSOCIATED WITH THE STIFF.NESS DISTRIBUTION 

The increase in structural weight associated with a given 
increase in stiffness can be estimated on the basis of assump- 
tions similar to those made in relating the stiffness and the 
strength. For the purpose of this analysis the various com- 
ponents of the wing structure, exclusive of the carry-through 
structure within the fuselage, are classified in two groups: 
one which contains tlie elements that take the bending and 
torsional loads due to the assumed loading and one which 
contains all other components. In the first group are 

( 1 ) The amount of top and bottom skin that is used in the 
estimation of the thicknesses required to withstand the 
bending and torsional loads, including stringers and spar 
flanges included in the equivalent skin 

( 2 ) Webs, including any web stiffeners 
In the second group are 

( 1 ) Skin, stiffeners, false spars, and so on, which are not 
considered in the estimation of the equivalent thicknesses 

(2) Ribs, bulkbeads, and posts designed to raise the buck- 
ling strength of the cover sheets 

(3) Control surfaces and their supports, attachments, 
and actuating mechanisms 

(4) The supports of internal stores 

This analysis is concerned only with the first group and, 
more specifically, with the relative increase in the weight of 
this group occasioned by an increase in stiffness of the main 
structure. Means of estimating the actual magnitude of the 
weights involved and of estimating the weights of some of the 
items in the second group as well are given in references 9 
and 10 . 

The weight per unit length of the structural elements of the 
first group can be written as 

Wj = 27)87)2i7jCt« (B2S) 

where 7 , is the density of the material of the primary struc- 
ture (or an equivalent density in the case of sandwich con- 
struction), 1 J 21 is the ratio of an equivalent perimeter p to 
the actual perimeter of the cell, and p is the sum of all the 
lengths which constitute the perimeter, each multiplied by 
the ratio of its equivalent thickness to the equivalent 
thickness i, of the upper cover sheet. 

In view of the assumption made concerning the combina- 
tion of bending and torsional stresses, the thickness re- 
quired in equation (B28) can be obtained from equation 
(B18) as 

f =-"^-. 

' vuch^ 
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SO that 


or 




TUa ft' 


W;_ Ijlr 

Ws~{hfhr)- 


(B29) 


Consequently, the total weight (for both wings) of the struc- 
tural elements of the first group can be estimated from the 
relation 


V2 r* 

, 72 Jo 

cos , 


IIL 


ih/hrf 


ds*- 


(B30) 


Equation (B30) serves to estimate relative changes in tlie 
weight of the first group of structural elements. For instance, 
with a given distribution of / and h. that weight is directly 
proportional to C and inversely proportional to A/'. vSimi- 
larly. given two different distributions of / and It with the 
same values at the roof, the ratio of the wciglit.s is equal to 
the ratio of the two values obtained by using the res])ective 
distributions of I and A in the integral of e(|uation (B30). 

Although the actual value of IF, is not relevant to tliis dis- 
cussion, it may be estimated by substituting the previously 
calculated stiffness distributions into e(|ualion (B30), and 
tbe result is given here as a matt ei- of general interest: 


where 


V3VoVsV-2i Tv- It IF , b';2 j. 

Us — 7 , f ir 

’liiMo’Jis I n it^ COS A 

Cr 


F = — - F r* ch"' 

" (Id-xy^ Jo (A/A.,)'-’"^ 


(B31) 


(B32) 


According to equation (B31), the sti-uclui-al weight is di- 
rectly proportional to the design gross weight, load factor, 


swept-span aspect ratio, span, and density of the material 
of the primary structure and inversely proportional to the 
allowable stress and the wing thickness ratio. The dejicnd- 
ence of the weight on the taper ratio (all other parameters, 
notably tbe aspect ratio and span, are the same) is illustrated 
in figure 14 by a plot of the function against taper ratio for 


several values of the parameter and for several ratios of 
the wing thickness ratios at the tip and at the root. 
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